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Аннотация 

Введение. Математическое моделирование эффективно при анализе промышленной безопасности на 
металлургических предприятиях, в частности для отслеживания проблем системы «человек — машина». Чтобы 
ввести в рассмотрение фактор времени, задействуют рекуррентные соотношения (в дискретной модели) и 
дифференциальные (в непрерывной). Однако необходимо также решить проблему привязки параметров модели 
к реальным условиям производственной среды и к человеческому фактору. Цель данного 
исследования — создание метода определения параметров имитационных математических моделей динамики 
психофизиологических показателей оператора, влияющих на его работу. 

Материалы и методы. Психофизиологическое состояние (ПФС) оператора оценивали по работоспособности, 
утомляемости и ошибаемости. Данные собрали по тесту цифровой корректурной пробы (ЦКП). На основании 
полученных результатов вычислили экспериментальные значения показателей ПФС оператора, которые привели 
к нормированной шкале [0, 1]. Эти показатели для конкретного респондента, математическую модель и 
разработанный алгоритм задействовали при определении числовых значений параметров модели. Для 
интерпретации показателей работоспособности, утомляемости и ошибаемости ввели шкалы с пятью градациями. 
Результаты исследования. Использование модифицированного авторами варианта математической модели 
показало значительное улучшение ее прогностических свойств. Из 10 участников наилучший результат оказался у 
респондента № 7, худший — у респондента № 8. В течение 1-го часа работы (с 9.00 до 10.00) их работоспособность 
выросла примерно одинаково, с 0,5—0,55 почти до 0,6. Затем показатель респондента № 7 активно увеличивался и до 
конца рабочего дня оставался существенно выше уровня «хороший». Показатель респондента № 8 падал и с 14.00 до 
15.00 оказался ниже среднего. Разницу во многом определили хронотипы операторов. Их хронофизиологические 
особенности сказались также на утомляемости и ошибаемости. Для разных участников экспериментов 
варьировалось качество модели. В одном случае оно оказалось отличным (средняя относительная ошибка <5 %), в 
трех случаях — хорошим (<10 %), в четырех — удовлетворительным (<15 96). 

Обсуждение и заключение. Предлагаемый подход дает возможность получить для каждого индивидуума 
динамические профили его психофизиологических характеристик, оценить их взаимосвязи и выполнить прогноз 
на основе модифицированной математической модели. Однако для расширения функциональных возможностей 
моделей в реальных условиях работы оператора металлургического оборудования требуется увеличить объем 
выборки, уменьшить шаг дискретного времени и выполнить исследования для различных условий работы с 
учетом технологических, климатических, экологических, психологических и прочих факторов. 


Ключевые слова: безопасность на металлургических предприятиях, система «человек — машина», 
психофизиологическое состояние оператора, хронотип оператора, работоспособность оператора 
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Abstract 

Introduction. Mathematical modelling is effective in the analysis of industrial safety at metallurgical plants, in 
particular for tracking problems of the man — machine system. To introduce the time factor, recurrence relations (in a 
discrete model) and differential relations (in a continuous model) are used. However, it is also necessary to solve the 
problem of linking the model parameters to the real conditions of the production environment and to the human factor. 
The aim of this study is to create a method for determining the parameters of simulation mathematical models of the 
dynamics of the operator’s psychophysiological indicators affecting the work. 

Materials and Methods. The operator’s psychophysiological state (PPS) was assessed by performance, fatigue levels, 
and error rate. The data were collected by the Digital Correction Task (DCT) test. Based on the obtained results, the 
experimental values of the operator’s PPS indicators, which were reduced to the normalized scale [0, 1], were 
calculated. These indicators for a particular respondent, the mathematical model and the developed algorithm were used 
to determine the numerical values of the model parameters. In order to interpret the indicators of performance, fatigue 
and error rate, we introduced scales with five gradations. 

Results. The use of the authors’ modified version of the mathematical model showed a significant improvement in its 
prognostic properties. Out of 10 participants the best result was shown by respondent no. 7, the worst result was shown 
by respondent no. 8. During the first working hour (from 9.00 to 10.00) their performance increased almost equally, 
from 0.5—0.55 to almost 0.6. Then the score of respondent no. 7 increased and remained well above the “good” level 
until the end of the day. The score of respondent no. 8 dropped and was below average from 14.00 to 15.00. The 
difference was largely determined by the operators’ chronotypes. Their chronophysiological characteristics also affected 
fatigue and error rate. The model’s quality varied for different participants in the experiments. In one case it was 
excellent (mean relative error <5%), in three cases it was good (<10%) and in four it was satisfactory (<15%). 
Discussion and Conclusion. The proposed approach allows us to obtain the dynamic profiles of psychophysiological 
characteristics for every individual, to assess their interrelationships and to perform a prediction on the basis of a 
modified mathematical model. However, in order to extend the functionality of the models to the real working 
conditions of the metallurgical plant operator, it is necessary to increase the sample size, reduce the discrete time step 
and conduct studies for different working conditions, considering technological, climatic, environmental, psychological 
and other factors. 


Keywords: safety at metallurgical enterprises, man — machine system, psychophysiological state of the operator, 
operator chronotype, operator performance 
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Введение. Согласно данным Федеральной службы по экологическому, технологическому и атомному 
надзору Российской Федерации, к началу 2022 года на металлургических и коксохимических производствах 
зарегистрировано 1280 опасных производственных объектов, 19 из которых относятся к первому классу 
опасности, 325 — ко второму, и 922 — к третьему!. Как пишет журнал «Промышленные страницы»”, в 
2022 году на российских металлургических предприятиях самыми частыми происшествиями были пожары и 
взрывы. Эти инциденты обусловлены техническими причинами и нарушениями правил промышленной 
безопасности. Например, в январе 2022 года трое сотрудников Новолипецкого металлургического комбината 
проводили плановый осмотр межцехового газопровода без обязательной газозащитной аппаратуры. Рабочие 
погибли в результате отравления токсичными парами и газом. 

В [1] выделены технические, санитарно-гигиенические, организационные и психофизиологические группы 
факторов аварийности и травматизма на металлургических и коксохимических предприятиях. Последние две 
группы тесно связаны с человеком — участником производственного процесса. С точки зрения 
производственной безопасности металлургических предприятий особая ответственность возлагается на 
операторов высокотехнологичных агрегатов. В инженерной психологии их принято называть 
«человек-оператор». В рамках представленной работы в этом значении будем употреблять слово «оператор». 
Иногда именно он предотвращает переход опасности в инцидент или аварию на производстве. 

ГОСТ 12.0.003-20153 выделяет в самостоятельный блок опасные и вредные производственные факторы 
психофизиологического воздействия на человека. Оператор испытывает нервно-психические перегрузки, 
связанные с напряженностью трудового процесса. В их числе: 

— умственное перенапряжение; 

— перенапряжение анализаторов; 

— монотонность труда; 

— неуверенность в действиях из-за недостатка образования и опыта. 

Эмоциональные перегрузки ведут к переутомлению, плохому самочувствию, стрессовому состоянию и T. п. [1]. 

Надежность человека как элемента сложной технической системы зависит от внутренних и внешних 
условий, которые меняются во времени. Есть еще одна переменная — сам человек. Работу выполняют люди с 
разными личностными качествами, здоровьем, опытом и пр. 

В [2] исследовались сведения о несчастных случаях (НС) в кузнечно-прессовом производстве (КПП) за 
1985, 1987 и 1989 гг. Отмечены некоторые особенности в распределении по количеству часов от начала работы 
до инцидента. Эти закономерности можно объяснить человеческим фактором, точнее суточным ритмом 
трудовой деятельности (рис. 1). 
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Рис. 1. Распределение НС в кузнечно-прессовом производстве по количеству часов от начала работы до НС [2] 


! Годовой отчет о деятельности Федеральной службы по экологическому, технологическому и атомному надзору в 2021 году. Москва: 2022. 
С. 407. URL: https://www.gosnadzor.ru/public/annual_reports/ogoBoi%200T4eT%203a%202021%20r.pdf (дата обращения: 22.11.2023). 

2 Опасная тенденция: аварии и ЧС на металлургических производствах в 2022-2023 годах. Промышленные страницы. URL: 
https://indpages.ru/safe/avareeee-na-myetallurgeechyeskeeh-proeezvodstvah/?ysclid=Inshtjh0z2208396380 (дата обращения: 22.11.2023). 

З ГОСТ 12.0.03-2015. Система стандартов безопасности труда. Опасные и вредные производственные факторы. Электронный фонд 
правовых и нормативно-технических документов. URL: https://docs.cntd.ru/document/1200136071 (дата обращения: 21.11.2023). 
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На графиках рис. 1 достаточно выражена синхронность пиков, а также минимальная аварийность B 
последние часы работы. Значит, время работы оператора — существенный фактор возникновения несчастного 
случая на производстве. 

В [3] экспериментально установлена необходимость учитывать хронотип оператора при определении 
профессионально важных качеств. 

Для предприятий металлургической отрасли методы математического моделирования эффективны при 
анализе ситуаций, которые: 

— связаны с промышленной безопасностью; 

— могут возникнуть в сложной технической системе «человек — машина»; 

— формируются и развиваются в различных условиях внутренней и внешней среды. 

Математические модели, описывающие психофизиологическое состояние оператора, разрабатывались в 
различных математических постановках. При этом задействовали инструментарий измерения первичных 
показателей. В.Г. Абашин выяснил, как использование биометрических технологий связано с 
психофизиологическим состоянием оператора автоматического рабочего места, его работоспособностью, 
снижением количества брака и аварий технологического процесса по вине человека [4]. Работоспособность 
оператора моделируется по клавиатурному почерку. С этой целью применяют нечеткие множества как основу 
адаптивных моделей — биометрической и мультибиометрической [5]. 

Авторы [6] используют в математической модели биометрические признаки (голос, клавиатурный почерк и 
характер работы с компьютерной мышью). Это позволяет судить о психофизиологическом состоянии 
оператора: норма, усталость, опьянение, возбуждение, расслабление (засыпание). Модели базируются на 
стратегии Байеса, а также нейросетевом подходе и позволяют оценить уровень ПФС оператора и 
спрогнозировать его способность к выполнению текущих операций. 

В [7] разработаны имитационные математические модели взаимосвязи факторов психофизиологического 
состояния оператора. Авторы задействовали рекуррентные соотношения в дискретной модели и системы 
дифференциальных уравнений в непрерывной модели, что позволяет ввести в рассмотрение фактор времени. В 
этом случае параметры модели определяются по: 

— результатам тестирования конкретного работника; 

— его способности к определенному виду деятельности; 

— итогам экспериментальных исследований рабочего места и функциональных действий оператора. 

Эта часть моделирования меньше других формализуется и в наибольшей степени влияет на возможности 
использования модели для корректировки реальных действий оператора. В настоящее время проблема не 
решена. Таким образом, сохраняют актуальность вопросы определения параметров математических моделей, 
которые могут корректно описать взаимосвязи различных факторов психофизиологического состояния 
оператора. Цель данной статьи — разработка методики определения параметров таких моделей. Научные 
изыскания в данной сфере будут способствовать выявлению показателей состояния оператора 
металлургического оборудования, которые влияют на его функциональные характеристики и могут стать 
причиной серьезных ошибок с негативными последствиями. 

Материалы и методы. Для оценки психофизиологического состояния оператора выбрали следующие 
показатели: 

— работоспособность (возможность целенаправленно выполнять работу определенного количества и 
качества за обозначенное время); 

— утомляемость (снижение работоспособности с нарушением координации движений, падением 
концентрации внимания и точности решений [8]); 

— ошибаемость (оценивается количеством ошибочных действий). 

Утомление определяется по субъективным и объективным признакам [9]. Субъективный уровень утомления 
устанавливали по шкале оценки усталости“ (fatigue assessment scale, РАЗ) [10]. Ее разработала группа 
нидерландских ученых под руководством X. Михильсен (H.J. Michielsen). Это 10 вопросов о ежедневном 
самочувствии респондента. 

В качестве источника объективных параметров в данной работе использовали тест цифровой корректурной 
пробы (ЦКП) По его результатам вычислялись экспериментальные значения показателей 
психофизиологического состояния оператора. Числовые значения параметров модели определяются по 
названным показателям для конкретного респондента в совокупности с используемой моделью. 


* Шкала оценки усталости: Fatigue Assessment Scale (FAS). URL: https://wasog.org/dynamic/media/78/documents/Questionairres/fas rus anon.htm] (дата 
обращения: 21.11.2023). 
5 Цифровая корректурная проба. URL: https://metodorf.ru/tests/korrekt/korrektchis.php (дата обращения: 18.08.2023). 
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Задача подбора параметров математической модели взаимосвязи различных психофизиологических 
факторов оператора решается в два этапа. 

1-й — экспериментальный. Участников тестируют, собирают данные об условиях работы оператора и его 
личностных качествах, которые могут влиять на исследуемые параметры трудовой деятельности. 

2-й — расчетный. Используется база данных, собранных на 1-м этапе. С учетом выбранной математической 
модели оцениваются ее параметры для каждого конкретного работника. При достаточном объеме 
статистических данных в широком диапазоне варьируемых параметров строятся усредненные модели, которые 
можно применить для прогнозирования изменения психофизиологических показателей работы оператора в 
различных условиях и с учетом временного фактора. 

В данном исследовании для электронного тестирования использовали программу, представленную 
интернет-ресурсом «Интерактивный портал — книга методик саморазвития и достижения успеха». 

При выполнении ЦКП испытуемый строка за строкой просматривал сгенерированный программой числовой 
массив и зачеркивал цифры, указанные в задании. На выполнение теста давали 3 минуты. Затем программа 
выдала описанные ниже результаты. 

1. Основные первичные показатели: 

— время выполнения теста f; 

— общее количество просмотренных цифр до последней выбранной цифры № 

— общее количество просмотренных строк С; 

— общее количество цифр, которые необходимо вычеркнуть и; 

— общее количество вычеркнутых цифр М; 

— число верных ответов 5; 

— число пропущенных цифр Р; 

— число ошибочно выбранных цифр О. 

2. Расчетные показатели, характеризующие: 

— скорость (производительность) внимания A; 

— точность работы 7 (в трех вариантах); 

— умственную продуктивность Е; 

— умственную работоспособность Ан; 

— концентрацию внимания К; 

— устойчивость концентрации внимания Ky; 

— объем зрительной информации V; 

— скорость переработки информации О. 

Для решения поставленной задачи отобрали показатели, на основе которых оценивались 
работоспособность X, утомляемость Уи ошибаемость Z оператора по формулам (1—4). 


1. Умственная работоспособность А, (по цифровой корректурной пробе”, единица измерения — знаков в секунду): 
A, =(N/t)x((M-(O+P))/n). (1) 


Для приведения показателя к нормализованной шкале [0, 1] в безразмерных единицах авторы предлагают 
формулу (2). Результатом будет i-Toe значение показателя работоспособности X (соответствует дискретному 
условному времени i): 


х = A, [тах { А, } > (2) 


где max {Áu} — максимально возможное значение умственной работоспособности по данному виду теста 
независимо от тестируемого работника. 

Рассмотрим стандартный тест ЦКП с 1600 цифрами. Допустим, респондент за 180 секунд полностью и 
правильно выполнит тест. В этом случае величина тах{А,} составит 8,889 знаков в секунду. Если же 
тестируемый выполнит тест менее чем за 180 с, значение X при расчете по формуле (2) может оказаться 
больше 1. Тогда для Х принимается значение 1. 

2. Утомляемость У оценивается по показателю К (концентрация внимания): 


Ү=(1-К /100), (3) 


где К = (М-0) п 100 — коэффициент, характеризующий концентрацию внимания, %. 


* Интерактивный портал — книга методик саморазвития и достижения успеха. URL: https://metodorf.ru/ (дата обращения: 18.10.2023). 
7 Цифровая корректурная проба. URL: https://metodorf.ru/tests/korrekt/korrektchis.php (дата обращения: 18.10.2023). 
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Данный показатель У принимает безразмерные значения в интервале [0, 1]. 
3. Ошибаемость Z предлагается оценить по формуле: 


2=(0+Р)/(М +P). (4) 

Данный показатель также принимает значения в интервале [0, 1]. 
В дальнейшем будем характеризовать X, У и 7 как унифицированные количественные индикаторы 
психофизиологического состояния оператора, связанные с его основной производственной деятельностью. Для 


трактовки показателей X, Уи Z разработана шкала, представленная в таблице 1. 


Таблица 1 
Интерпретация значений показателей в качественной форме 
Показатель Интервалы значений показателей / качественная характеристика 
[0-0,20] [0,20—0,40] [0,40—0,60] [0,60—0,80] [0,80—1,00] 
Работоспособность X 
низкая ниже средней средняя хорошая высокая 
[0—0,20] [0,20—0,40] [0,40—0,60] [0,60—0,80] [0,80—1,00] 
Утомляемость Y 
низкая ниже средней средняя выше средней высокая 
[0—0,01] [0,01—0,05] [0,05—0,10] [0,10—0,20] [0,20—1,00] 
Ошибаемость Z 
незначительная заметная существенная значительная высокая 


Показатели X, Уи 7 — безразмерные, унифицированные со значениями в интервале [0, 1]. 

На втором этапе для подбора параметров математической модели есть два способа. 

1. При помощи модуля «Поиск решения» программы Excel создается сценарий решения системы для 
некоторых начальных параметров. В его основе — система рекуррентных уравнений математической модели в 
дискретной форме. Затем «Поиск решения» выдает оптимальные значения параметров модели для целевой 
функции, равной сумме квадратов отклонений расчетных значений показателей X, Y и Z от экспериментальных. 

2.В программе Anylogic?, которая позволяет в автоматическом режиме подбирать параметры 
рассматриваемых математических моделей в дискретной форме рекуррентных соотношений и в виде системы 
дифференциальных уравнений.  Anylogic целесообразно использовать после анализа результатов 
предварительной оценки параметров модели с помощью модуля «Поиск решения» Excel. 

Для апробации предлагаемой методики отобрали 10 человек — студентов и сотрудников Донбасского 
государственного технического университета. Они работали как операторы на учебных тренажерах с 
автоматическим и полуавтоматическим управлением. В начале и в конце рабочего дня оценивался 
субъективный уровень утомления по шкале FAS. В начале каждого часа, с 9.00 до 16.00, выполнялись тесты 
ЦКП. Данные обрабатывали по формулам (1-4) для каждого респондента. Результаты тестирования и их 
унификацию (X, Уи Z) использовали на втором этапе для определения параметров математической модели. 


В данной работе применили модель в рекуррентной форме [7], которая описывается системой уравнений: 


Ха = № + а ^ a ^i) byx;y; — hx;z;, 
dix; + с, К 
Ji Yi 
ы = Y; + ——— (1-2) + x,y, 5 
Ян = У "dy +6, ( my оу (5) 
7. Lon 
Zj41 = 2; +a +. (1-—+) + by y,z;. 
+1 3 daz, +e, ( E зу 


Решением системы (5) являются три сопряженных временных ряда длиной т индикаторов 
X = (x, i=0, 1, 2, ..., т}, Y= {у,, 1=0, 1,2,....т}, Z= {z,i=0, 1, 2,..., m}. Индекс: — переменная, 
характеризующая дискретное время в системе. Параметры aj bj К, h, 4, с; ( = 1, 2, 3) определяются в 


результате решения оптимизационной задачи по исходным данным тестирования. 


$ AnyLogic: имитационное моделирование для бизнеса. URL: https://www.anylogic.ru/ (дата обращения: 12.11.2023). 
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Математическая модель оптимизационной задачи описывается системой ограничений (6) и целевой 
функцией (7). 
o, (aj, b, kde; ) = (xp -xf P SA, 
в -») < А, 
8,;(a;,b;,k;,h,d,,c;)=(2? -z P <A, 
a,,b;,k;,h,d;,c; 20, (6) 
xP 20, у? 2 0,z? > 0, 


xP < Ту? <1,2Р <1, 


i-l..,m, 
j=1,2,3. 
V(a, bj kj hdc) - S B, € 9 By + 9 5, > min. (7) 
i=l i=l i=l 


Здесь индексом ф помечены фактические значения индикаторов X, Y и Z, найденные B результате обработки 
данных тестирования оператора, р — расчетные значения X, Y и 7, определяемые из решения системы (5) и 
зависящие OT искомых в данной задаче оптимизации переменных ау, by, kj, h, dj, cj (ў = 1, 2, 3). Ha варьируемые 
переменные а, bj kj, h, dj, с; ограничения задаются по максимально допустимому значению А квадратов 
отклонений (Öxi, бу, бы) расчетных значений X, Уи Z от фактических для всех т значений рассматриваемых 
временных рядов. Целевая функция (а, bj, kj, h, dj, cj) равна сумме квадратов отклонений (Öxi, бук, Өш). 

Для решения оптимизационной задачи применялся метод обобщенного приведенного градиента (ОПГ). 

Результаты исследования. 1-й этап — экспериментальный. Результаты тестирования участников 
эксперимента (респонденты №1-№10), преобразованные к унифицированным показателям X, Уи Z, представлены 
на диаграммах размаха (рис. 2-4). В виде точки показан центр распределения (медиана), прямоугольник 
обозначает границы варьирования (квартили 25 % — 75 %), «усики» — нижнюю и верхнюю границу значений 


показателя (min — тах). Диаграммы выполнены с помощью программного комплекса «Statistica». 
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Рис. 2. Показатели работоспособности операторов 
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Рис. 3. Показатели утомляемости операторов 
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Рис. 4. Показатели ошибаемости операторов 


Анализ рис. 2—4 позволил выделить респондентов с лучшими и худшими оценками их деятельности 
(таблица 2). Учитывались уровни показателей X, Уи Z, соотнесенные с градациями таблицы 1, а также квартили 


и интервалы варьирования показателей. 
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Таблица 2 
Ранжирование респондентов по результатам тестирования 
E Варианты (Ne респондента) 
Критерий оценки 
наихудшие наилучшие 

Работоспособность Х (чем больше, тем лучше) 8,1,4,6 7,9,3 
Утомляемость Y (чем больше, тем хуже) 8,1,4,6 7,10 
Ошибаемость Z (чем больше, тем хуже) 4,8,1 9, 10, 7 
Интегрально (с учетом Х, У, Z) 8,4, 1 7, 10, 9 


Графики на рис. 5 отражают динамику работоспособности для респондентов №7 (пунктирная линия, 
лучший результат) и № 8 (сплошная линия, худший результат). 
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Рис. 5. Динамика работоспособности операторов 


Для объяснения достаточно большой разницы данных респондентов 7 и 8 учитывали их возраст, пол, род 
занятий, самочувствие и др. Наиболее существенным оказался фактор хронотипа. Классификация, принятая в 
1970 году, подразумевает три категории людей с отличительными поведенческими признаками и генетической 
разницей биоритмов. Это так называемые «жаворонки», «совы» и «голуби» [11]. Респондент № 8по 
хронотипу — «жаворонок», №7 — «сова». В результате работоспособность № 8 снижается к концу дня, 
а №7 — повышается (рис. 5). Хронофизиологические особенности операторов проявляются также в динамике 


их утомляемости и ошибаемости. 


2-й этап — расчетный. Получены параметры модели (5), рассчитанные для каждого участника 
эксперимента. С этой целью решили оптимизационную задачу (6)-(7). Диапазоны изменения параметров 


приведены в таблице 3. 


Таблица 3 


Диапазоны параметров математической модели (5) для группы участников эксперимента 


Показатель Диапазоны параметров Показатель Диапазоны параметров 

a1 0,0004—0,0300 bs 0,0001—10,1400 
bi 0—0,9700 Кз 0,0020—0,0700 
ky 0,5900—1,0000 Сі 0,0010-0,2000 
h 0—1,8100 с? 0,0001—0,0010 
a2 0—0,0010 C3 0,0100—0,3300 
b» 0,0070—0,2600 di = 0 

ky 0,4200—1,0000 dy =0 

аз 0,0005—0,0600 ds 0,0400—0,9800 


B таблице 4 представлены результаты оценки качества математических моделей динамики изменения 
показателей X, Уи Z, оцениваемых по среднеквадратической ошибке модели и средней относительной ошибке 


для всех испытуемых. 
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Таблица 4 
Оценка качества математических моделей 


№ респондента 1 2 3 4 5 6 7 8 9 10 


Средняя относительная ошибка =, % | 14,27 | 4,54 | 20,93 | 12,80 | 7,28 | 12,00 | 9,66 | 15,03 | 12,23 | 9,70 


Среднеквадратическая ошибка, ед. | 0,182 | 0,107 | 0,253 | 0,267 | 0,158 | 0,247 | 0,271 | 0,256 | 0,192 | 0,166 


Для разных респондентов показатели качества модели варьируются в широком диапазоне. Отличное 
качество (= < 5 %) получено в одном случае, хорошее (= < 10 %) — в трех случаях, удовлетворительное 
(= < 15%) — в четырех. Для двух случаев (№3 и №8) не удалось удовлетворительно решить задачу 
оптимального подбора параметров математической модели. Это говорит о том, что не все влияющие факторы 
учтены либо модель (5) в некоторых случаях не работает. 

Согласно результатам анализа алгоритма решения поставленной задачи, ряд упрощений улучшит 
сходимость результатов. Можно, например: 

— уменьшить на 3 ед. число параметров системы (5), сократив числители и знаменатели вторых слагаемых 
правых частей системы на C1, C2, сз соответственно; 

— удалить из системы (6) ее первое, второе и третье ограничения, переведя их в статус наблюдаемых ограничений. 

В итоге алгоритм метода ОПГ будет работать лучше. 

На рис. 6 и 7 приведены сравнительные диаграммы оценки качества построенных математических моделей 
в двух вариантах: 

— 1 — исходная модель (5); 

— 2 — модифицированная модель. 

Как видно из диаграмм, решение задачи подбора параметров математической модели динамики показателей 
ПФС оператора значительно улучшено при использовании модели 2 для всех случаев, кроме № 8. 
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Рис. 6. Сравнение математических моделей по средней квадратической ошибке 
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Рис. 7. Сравнение математических моделей по средней относительной ошибке 
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Обсуждение и заключение. Предложенный подход позволяет получить и прогнозировать для каждого 
индивидуума динамические профили его психофизиологических характеристик, составить математические 
модели взаимосвязей. Целесообразно увеличить объем выборки, уменьшить шаг дискретного времени и 
выполнить исследования для различных условий работы, включающие технологические, климатические, 
экологические, психологические и прочие факторы, которые могут влиять на эффективность работника. Это 
следует сделать для расширения функциональных возможностей моделей, их применения в реальных условиях 
работы оператора металлургического оборудования. 

Обнаруженный эффект влияния хронотипа человека на его ПФС требует более основательных 
исследований. В статье [12], посвященной поиску компонентов циркадных часов у людей, на больших 
статистических выборках показано, что хронотип зависит от многих факторов (пол, возраст, график работы и 
др.). Все это необходимо диагностировать и учитывать в математических моделях оценки ПФС операторов. 

Интересно также дополнить модель компонентами культуры производственной безопасности [13], которые 
можно привести к количественному виду и использовать как поправочные коэффициенты. 

В дальнейшем предполагается усовершенствовать модель дополнительными переменными. Для этого нужно 
зафиксировать психофизиологические показатели и определять местонахождения оператора в режиме 
реального времени. 

Представленную в работе методику можно взять за основу при решении описанных ниже задач. 

— Составление и анализ динамического профиля работника при приеме на работу в качестве оператора 
машины, агрегата или устройства, где важны такие характеристики ПФС, как работоспособность, утомляемость 
и ошибаемость. Эта задача решается для конкретного человека в определенных условиях производства, что 
позволяет подобрать оптимальный режим работы и отдыха оператора, способствует сохранению здоровья 
специалиста и повышению уровня промышленной безопасности предприятия. 

— Формирование и ведение специализированных баз статистических данных, включающих характеристики 
рабочих мест и сотрудников, в том числе их психофизиологическое состояние. На основе собранного 
статистического материала и систем мониторинга ПФС оператора в режиме реального времени можно строить 
прогнозные модели для предупреждения нештатных и аварийных ситуаций на металлургических предприятиях. 
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Аннотация 

Введение. Нормативными документами по пожарной безопасности при определении численности и 
технической оснащенности подразделений пожарной охраны в явном виде не учитываются природно- 
климатические и географические особенности субъектов Российской Федерации. При этом в научной 
литературе встречаются сведения о влиянии отдельных природно-климатических факторов на эффективность 
действий пожарной охраны при тушении пожаров. Целью настоящего исследования явилось определение 
влияния системы природно-климатических и географических характеристик субъектов Российской Федерации 
на показатели оперативного реагирования территориальных подразделений пожарной охраны. Полученные 
результаты рекомендованы для дальнейшего использования при нормировании численности и технической 
оснащенности подразделений пожарной охраны. 

Методы и материалы. Проведен факторный анализ статистических данных за 2020—2022 гг. Статистические 
данные получены из федеральной государственной информационной системы «Федеральный банк данных 
«Пожары», с сайта Федеральной службы государственной статистики Российской Федерации и из других 
источников. Для анализа отобраны 10 показателей, характеризующих природно-климатические и 
географические особенности субъектов Российской Федерации, 2 показателя оперативного реагирования 
территориальных подразделений пожарной охраны и 4 показателя обстановки с пожарами. 

Результаты исследования. Выделены пять значимых факторов, изменение которых объясняет изменение 
наблюдаемых показателей. Первый фактор характеризует связь климатических условий с показателями 
обстановки с пожарами. Второй фактор связывает показатели оперативного реагирования подразделений 
пожарной охраны с особенностями рельефа субъектов Российской Федерации. Третий фактор описывает 
взаимосвязь показателей обстановки с пожарами и показателей оперативного реагирования с плотностью 
населения и лесистостью территории. Остальные факторы в показатели обстановки с пожарами и оперативного 
реагирования существенного вклада не вносят. 

Обсуждение и заключение. С помощью математического анализа и с использованием факторного 
моделирования авторами исследована взаимозависимость природно-климатических и географических 
особенностей субъектов Российской Федерации, показателей обстановки с пожарами и показателей 
оперативного реагирования подразделений пожарной охраны. Определены наиболее значимые факторы, 
влияющие на эти показатели. Среди них средняя температура воздуха, лесистость территории, наличие горных 
массивов, плотность населения. Данные показатели следует учитывать при определении численности и 
технической оснащенности подразделений пожарной охраны для повышения эффективности их 
функционирования. 


Ключевые слова: пожарная охрана, факторный анализ, субъект Российской Федерации, гибель, травматизм, 
время прибытия, природно-климатические и географические показатели 
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Abstract 

Introduction. When determining the number and technical equipment of fire departments, regulatory documents on fire 
safety do not take into account natural, climatic and geographical features of the subjects of the Russian Federation in 
an explicit form. At the same time, there is some information in scientific literature about the influence of certain 
natural and climatic factors on the effectiveness of fire protection actions in extinguishing fires. The aim of this study 
was to determine the influence of the system of natural, climatic and geographical characteristics of the subjects of the 
Russian Federation on the rapid response indicators of territorial fire departments. The results obtained are 
recommended for further use in rationing the number and technical equipment of fire departments. 

Materials and Methods. A factor analysis of statistical data for 2020—2022 was carried out. The statistical data were 
obtained from the Federal State Information System "Federal Database "Fires", from the website of the Federal State 
Statistics Service of the Russian Federation and from other sources. For the analysis, the authors selected ten indicators 
characterizing natural, climatic and geographical features of the subjects of the Russian Federation, two indicators of 
the rapid response of territorial fire departments and four indicators of the fire situation. 

Results. Five significant factors were identified, the change of which explained the change in the observed indicators. 
The first factor characterized the relationship of climatic conditions with fire situation indicators. The second factor 
connected the indicators of the rapid response of fire departments with the terrain features of the subjects of the Russian 
Federation. The third factor described the relationship between fire situation indicators and rapid response indicators 
with population density and forest cover of the territory. Other factors did not significantly contribute to the indicators 
of fire situation and rapid response. 

Discussion and Conclusion. By means of mathematical analysis and factor modeling, the authors investigated the 
interdependence of natural, climatic and geographical features of the subjects of the Russian Federation, fire situation 
indicators and indicators of rapid response of fire departments. The most significant factors influencing these indicators 
were identified. They included the average air temperature, the area covered by forest, the presence of mountain ranges, 
and population density. These indicators should be taken into account when determining the number and technical 
equipment of fire departments to increase the efficiency of their functioning. 


Keywords: fire protection, factor analysis, subject of the Russian Federation, death, injury, arrival time, climatic and 
geographical indicators 
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Введение. При определении численности и технической оснащенности территориальных подразделений 
пожарной охраны проводится обследование территории субъекта Российской Федерации, которое включает в 
себя анализ состояния существующей системы обеспечения пожарной безопасности субъекта Российской 
Федерации, оперативно-тактических характеристик (особенностей) рассматриваемой территории, частоты 
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возникновения пожаров, рисков гибели и травмирования на них людей'. Среди факторов, влияющих на 
состояние пожарной безопасности в субъектах Российской Федерации, большое значение имеют природно- 
климатические и географические характеристики. В этой связи в настоящей работе с использованием 
факторного анализа было проведено исследование взаимосвязи природно-климатических и географических 
характеристик субъектов Российской Федерации, показателей оперативного реагирования территориальных 
подразделений пожарной охраны и показателей обстановки с пожарами. Факторный анализ является одним из 
наиболее широко используемых методов, направленных на выделение из множества разнородных показателей 
небольшого количества латентных факторов, вариация которых может объяснить изменение большей части 
наблюдаемых показателей [1—5]. Это позволяет построить математическую модель, в которой исследуемые 
факторы имеют простую и наглядную структуру [6]. Ранее авторами с использованием факторного анализа 
проводилось изучение оперативного реагирования подразделений пожарной охраны [7], социальных 
последствий пожаров [8], обстановки с пожарами в сельских поселениях [9], готовности подразделений 
пожарной охраны [10]. 

Материалы и методы. Проводимое исследование включало в себя три этапа. На первом этапе была 
сформирована матрица исходных показателей. Для этого в каждом субъекте Российской Федерации были 
определены значения 16 показателей, разбитых на три группы. 

В первую группу вошли данные, характеризующие обстановку с пожарами (обозначение показателя указано 
в скобках): 

— количество пожаров, приходящихся на | тыс. человек населения, ед. (71); 

— средний прямой материальный ущерб от одного пожара, руб. (Y); 

— количество погибших при пожарах, приходящихся на 100 тыс. человек населения, чел. (Үз); 

— количество травмированных при пожарах, приходящихся на 100 тыс. человек населения, чел. (Y4). 

Во второй группе отражались показатели оперативного реагирования подразделений пожарной охраны: 

— среднее время, за которое прибывает первое подразделение пожарной охраны, мин. (71); 

— среднее время, затрачиваемое на тушение пожара, мин. (Zo). 

Третью группу составляли показатели, характеризующие природно-климатические, географические и 
демографические особенности субъектов Российской Федерации: 

— плотность населения, чел/км? (X1); 

— доля водной поверхности от площади территории субъекта, % (X2); 

— длина морской береговой линии от общей протяженности границы субъекта, % (X3); 

— доля лесов от площади территории субъекта, % (X4); 

— среднегодовое количество осадков, MM в год (X5); 

— средняя температура июля, °С (X6); 

— средняя температура января, °С (X7); 

— сейсмическая опасность (наличие населенных пунктов с указанной сейсмической интенсивностью для 
степени сейсмической опасности), °С (Хз); 

— доля территории, занятой горами, от общей площади территории субъекта, % (Xo); 

— количество солнечных дней в году, (X10). 

Показатели обстановки с пожарами и показатели оперативного реагирования подразделений пожарной охраны 
получены из электронной базы данных пожаров, которая ведется с использованием автоматизированной 
аналитической системы поддержки и управления контрольно-надзорными органами МЧС России?. 

Сведения о плотности населения и площади лесов в субъектах Российской Федерации взяты из данных 
Федеральной службы государственной статистики Российской Федерации?. Характеристики климата субъектов 
Российской Федерации получены на справочно-информационном портале «Погода и климат»“. При 
определении сейсмической опасности в субъектах Российской Федерации использованы данные 
СП 14.13330.20145. 


! Об утверждении методик расчета численности и технической оснащенности подразделений пожарной охраны. Приказ МЧС России or 
15.10.2021 № 700. URL: https://base.garant.ru/403 136953/ (дата обращения: 18.12.2023) 

2 06 утверждении Регламента работы в информационной системе «Автоматизированная аналитическая система поддержки и 
‘управления контрольно-надзорными органами МЧС России». Приказ МЧС России № 954 от 04.10.2022. URL: https://fireman.club/normative- 
documents/prikaz-mchs-rossii-954-ot-04- 1 0-2022-ob-utverzhdenii-reglamenta-raboty-v-informaczionnoj-sisteme/ (дата обращения: 11.12.2023) 

3 Регионы России. Социально-экономические показатели 2022. Федеральная служба государственной статистики Российской Федерации. 
URL: https://rosstat.gov.ru/folder/210/document/13204 (дата обращения: 08.12.2023) 

^ Справочно-информационный портал «Погода и климат». URL: http://www.pogodaiklimat.ru/ (дата обращения: 08.12.2023) 

5 Свод правил. Строительство в сейсмических районах. СП 14.13330.2018 от 25.11.2023. URL: https://docs.cntd.ru/document/550565571 (дата 
обращения: 08.12.2023) 
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Отобранные статистические показатели характеризовались большим разбросом, что было обусловлено 
различными размерами субъектов Российской Федерации (площадь территорий и численность населения), а 
также различиями их природно-климатических и географических условий. В связи с этим для приведения 
статистических данных к единообразной форме была проведена нормировка ряда показателей, что позволило 
преобразовать множество натуральных показателей в синтетические. Нормировка показателей осуществлялась 
на площадь территории субъектов (плотность населения, доля лесов и водной поверхности от площади 
территории субъекта) и на численность населения (количество пожаров, число погибших и травмированных 
при пожарах). 

Дальнейшее моделирование проводилось с использованием полученной матрицы синтетических 
показателей. Значимые факторы выделялись с использованием трех наиболее распространенных методов: 
центроидного, главных компонент и максимального правдоподобия. Все три метода дают близкие результаты, 
которые в пределах статистических погрешностей согласуются между собой. В дальнейшем для проведения 
факторного анализа был использован метод главных компонент. 

На втором этапе с помощью метода главных компонент были получены значимые факторы. Для облегчения 
их интерпретации выполнено вращение полученных факторов с использованием следующих методов: 

— варимакс (данный метод путем уменьшения числа переменных для каждого фактора позволяет выполнить 
лучшее разделение факторов); 

— квартимакс (данный метод путем уменьшения числа факторов, относящихся к каждой переменной, дает 
возможность выделить генеральный фактор и упростить интерпретацию); 

— биквартимакс; 

— эквимакс. 

Перед проведением процедуры вращения факторные нагрузки были нормализованы с использованием 
метода Кайзера. Это позволило исключить влияние переменных с большой схожестью на итоговый результат. 
В последующих исследованиях использовался метод вращения варимакс, поскольку анализ полученных 
результатов показал, что факторные нагрузки несущественно зависят от метода вращения. 

В результате факторного анализа было получено пять значимых факторов. Доля общей дисперсии, 
объясняемой данными факторами, распределена следующим образом: первый фактор — 22,1 %, второй 
фактор — 15,4 %, третий — 12,1 %, четвертый — 9,7 9o, пятый фактор — 7,7%. Данные пять значимых 
факторов объясняют 67,0 % общей дисперсии. 

В таблице 1 представлены значения факторных нагрузок для каждого из пяти факторов (значимые 
показатели выделены полужирным шрифтом). 


Таблица 1 
Матрица факторных нагрузок 
Переменная Фактор 
(показатель) 1 2 3 4 5 

Yı 0,019 -0,179 0,681 0,212 -0,308 

У 0,692 0,028 -0,241 0,027 0,123 

Y; 0,692 0,164 0,484 -0,123 0,021 

Y4 0,821 -0,095 0,229 0,046 -0,121 

21 -0,362 0,587 0,417 0,091 0,306 

% 0,009 0,218 0,509 -0,343 0,265 

X -0,093 0,069 -0,587 -0,086 -0,063 

X 0,018 0,234 0,056 0,070 -0,836 

X 0,000 -0,373 -0,155 -0,422 -0,532 

№ 0,236 -0,305 0,536 -0,357 0,002 

№ -0,334 0,026 0,008 -0,758 0,016 

Xs -0,636 0,032 -0,062 0,630 0,200 

X; -0,646 0,357 -0,343 -0,144 0,040 

Xe -0,073 -0,833 0,097 0,174 0,133 

Xo 0,058 -0,841 0,180 -0,019 0,051 

Xo -0,437 -0,322 0,191 0,630 -0,081 
Ha третьем этапе была решена творческая задача, которая выходит за пределы формального 
метода — проведена содержательная интерпретация полученных факторов с использованием предметных 

терминов. 
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Результаты исследования. Установлено, что в первый значимый фактор наибольший вклад вносят 
следующие переменные: 

— средний прямой материальный ущерб от одного пожара, руб. (Y2); 

— количество погибших при пожарах, приходящихся на 100 тыс. человек населения, чел. (Үз); 

— количество травмированных, приходящихся на 100 тыс. человек населения, чел. (Y4); 

— средняя температура июля, °С (X6); 

— средняя температура января, °С (X7); 

— количество солнечных дней в году (Х\о). 

Среди показателей обстановки с пожарами в первый фактор наибольший вклад вносит показатель 
«количество травмированных при пожарах, приходящихся на 100 тыс. человек населения». Данный фактор 
характеризует связь климатических условий с показателями обстановки с пожарами. В тех субъектах 
Российской Федерации, где более суровые климатические условия (низкие температуры, малое количество 
солнечных дней), отмечаются более тяжелые последствия пожаров (больше погибших и травмированных, выше 
материальный ущерб). 

Данная взаимосвязь демонстрируется на рис. | и 2, где показаны зависимости количества пострадавших 
(погибших и травмированных) при пожарах, приходящихся на 100 тыс. человек населения, от средней 
температуры января. Видно, что с уменьшением средней температуры увеличивается количество пострадавших 
при пожарах. 
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Puc. 1. Зависимость количества погибших при пожарах, приходящихся Ha 100 тыс. человек населения, OT средней 
температуры января в субъектах Российской Федерации. 
Прямая — результат аппроксимации линейной функцией у = -0,166х + 4,338 
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Puc. 2. Зависимость количества травмированных при пожарах, приходящихся 
на 100 тыс. человек населения, от средней температуры января в субъектах 
Российской Федерации. Прямая — результат аппроксимации линейной функцией 
у = -0,248х + 3,565 
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Для второго значимого фактора наиболее существенными являются следующие переменные: 

— среднее время, за которое прибывает первое подразделение пожарной охраны, мин. (71); 

— сейсмическая опасность (наличие населенных пунктов с указанной сейсмической интенсивностью для 
степени сейсмической опасности С) (Xs); 

— доля территории, занятой горами, от общей площади территории субъекта, % (Xo). 

Второй фактор характеризует взаимосвязь показателей оперативного реагирования подразделений 
пожарной охраны с географическими особенностями субъектов Российской Федерации. Наблюдается важная 
зависимость — чем большую территорию субъекта занимают горы, тем меньше среднее время, за которое 
первое подразделение пожарной охраны прибывает к месту вызова. 

Распределение факторных нагрузок для трех групп показателей в плоскости факторов | и 2 показано на рис. 3. 
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Рис. 3. Значения факторных нагрузок в плоскости факторов | и 2. Зеленым цветом обозначены показатели первой группы, 
красным — второй группы, синим — третьей группы 


В третий значимый фактор наибольший вклад вносят следующие переменные: 

— количество пожаров, приходящихся на | тыс. человек населения, ед. (Yi); 

— количество погибших при пожарах, приходящихся на 100 тыс. человек населения, чел. (Уз); 

— среднее время, за которое прибывает первое подразделение пожарной охраны, мин. (71); 

— среднее время, затрачиваемое на тушение пожара, мин. (72); 

— плотность населения, чел/км? (X1); 

— доля лесов от площади территории субъекта, % (X4). 

Третий фактор описывает взаимосвязь показателей обстановки с пожарами и показателей оперативного 
реагирования с плотностью населения и лесистостью территории. Как показали ранее проведенные 
исследования, с ростом плотности уменьшается средняя площадь области обслуживания одного подразделения 
пожарной охраны, результатом чего является снижение времени прибытия первого подразделения пожарной 
охраны. С другой стороны, чем ниже плотность, тем болыше площадь территории, на которой проживает 
одинаковое количество людей, соответственно, на большей площади происходит большее количество пожаров. 
Интересно также отметить, что с увеличением плотности населения снижается количество пострадавших при 
пожарах, приходящихся на 100 тыс. населения. 

На обстановку с пожарами также влияет площадь лесов. С увеличение доли территории субъекта 
Российской Федерации, занятой лесами, увеличивается количество пожаров, приходящихся в расчете на 
1 тыс. человек. 

Получена зависимость количества пожаров, приходящихся на | тыс. человек, от плотности населения. Как 
видно на рис. 4, количество пожаров уменьшается с ростом плотности населения. Данная зависимость 
описывается экспоненциальной функцией. 
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На рис. 5 представлено распределение факторных нагрузок для трех групп показателей в плоскости 
факторов | и 3. 
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Для четвертого значимого фактора наиболее важным являются следующие переменные: 

— среднегодовое количество осадков, мм B год (Xs); 

— средняя температура июля, C? (Xe); 

— количество солнечных дней в году (X10). 

Данный фактор связан с климатическими характеристиками субъектов Российской Федерации. Показатели 
первой и второй групп значительного вклада в данных фактор не вносят. Следует отметить, что взаимосвязь 
климатических характеристик субъектов с показателями обстановки с пожарами хорошо прослеживается из 
анализа факторных нагрузок для первого фактора (см. выше). 
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Для пятого значимого фактора определяющими являются следующие переменные: 

— доля площади водной поверхности от площади территории субъекта, % (X2); 

— длина морской береговой линии от общей протяженности границы субъекта, % (X3). 

Данные гидрографические характеристики не оказывают существенного влияния на обстановку с пожарами и 
на показатели оперативного реагирования подразделений пожарной охраны в субъектах Российской Федерации. 

Обсуждения и заключение. Впервые с использованием факторного анализа проведено исследование 
влияния природно-климатических и географических особенностей субъектов Российской Федерации на 
уровень пожарной безопасности. Проведенный анализ показал, что наиболее существенно на обстановку с 
пожарами и деятельность пожарной охраны в субъектах Российской Федерации влияют следующие природно- 
климатические, географические и демографические показатели: средняя температура, лесистость территории, 
наличие горных массивов, плотность населения. Таким образом, при определении численности и технической 
оснащенности подразделений пожарной охраны необходимо учитывать данные особенности субъектов 
Российской Федерации. 
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Аннотация 

Введение. Внутрилабораторные сличительные испытания являются важным и актуальным мероприятием для 
обеспечения качества и достоверности результатов исследований в лабораториях. Они позволяют оценить 
точность и воспроизводимость применяемых в лаборатории методов, а также выявить возможные источники 
ошибок и несоответствий в ее работе. Результаты внутрилабораторных сличительных испытаний 
предоставляются экспертам при подтверждении компетентности в рамках аккредитации. Как правило, 
сличительные испытания проводят в лабораторных условиях в привычной и спокойной для испытателей 
обстановке. Однако лаборатории, проводящие исследования в рамках специальной оценки условий 
труда (СОУТ), вынуждены проводить выездные сличительные испытания на реальных объектах, где за 
испытаниями наблюдают заказчики, которые невольно отвлекают работников лаборатории, что напрямую 
влияет на качество измерений. Целью данной работы является оценка качества проведения выездных 
внутрилабораторных сличительных испытаний на примере определения запыленности рабочего места 
каменщика на строительной площадке и определение минимально необходимого и достаточного для этого 
количества измерений. 

Материалы и методы. Для определения запыленности рабочего места каменщика использовался весовой 
метод, который заключается в сборе пыли на фильтры с последующим взвешиванием и определением 
концентрации. Оценка качества проведения внутрилабораторных испытаний производилась согласно 
ГОСТРИСО 5725-1-2002 «Точность (правильность и прецизионность) методов и результатов измерений. Часть 
1. Основные положения и определения». 

Результаты исследования. На рабочем месте каменщика выявлено превышение разовой предельно 
допустимой концентрации пыли в 1,6 раза. Средняя концентрация пыли на исследуемом рабочем месте 
составила: Km = 9,57+0,81 мг/м?, сходимость полученных результатов, r = 8,68 %, относительная погрешность, 
5 = 8,50 %. Выявлено, что максимально допустимая разница результатов двух испытаний составляет 0,84 мг/м}. 
Разница прямых измерений массы двух образцов должна быть не более 0,1 мг. 

Обсуждение и заключение. Полученные результаты показали возможность проведения минимального 
количества измерений, что по условиям воспроизводимости испытания признаются удовлетворительными и 
могут предоставляться экспертам для подтверждения компетентности лаборатории. В рамках СОУТ 


работникам исследуемого объекта рекомендуется использовать средства индивидуальной защиты органов 
дыхания, предохраняющие от высокодисперсной ПЫЛИ. 


Ключевые слова: сходимость, внутрилабораторный контроль, запыленность, строительная площадка, рабочее 
место, специальная оценка условий труда (СОУТ) 
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Abstract 

Introduction. Intra-laboratory comparison tests play an important role in ensuring the quality and reliability of research 
outcomes in laboratories. These tests allow researchers to evaluate the accuracy and reproducibility of the methods they 
use in their work, as well as to identify potential sources of error and inconsistency. The results of these tests are shared 
with experts to confirm competence within the accreditation. Typically, comparison tests are carried out in laboratory 
conditions in a familiar and calm environment for the testers. However, when laboratories conduct research as part of a 
special assessment of working conditions (SAWC) they are required to conduct on-site comparisons at real-world 
facilities, where customers may unwittingly disrupt the process and directly affecting the quality of the measurements. 
The aim of this study is to evaluate the quality of on-site intra-laboratory comparison tests using the example of 
determining the dust content in a bricklayer's work environment on a construction site, and to determine the minimum 
number of measurements necessary and sufficient for this purpose. 

Materials and Methods. To determine the dustiness of the bricklayer's workplace, a weighing method was used. This 
involved collecting dust on filters and then weighing them to determine the concentration. The quality assessment of 
intra-laboratory tests was conducted in accordance with GOST R ISO 5725-1-2002 "Accuracy (trueness and precision) 
of measurement methods and results. Part 1. General principles and definitions". 

Results. At the bricklayer's workplace, a 1.6-fold excess of the one-time maximum permissible dust concentration was 


detected. The average dust concentration at the workplace under study was: Кър = 9.57 + 0.81 mg/m’, the convergence 
of the results obtained was r= 8.68 %, the relative error б = 8.50 %. It was revealed that the maximum allowable 
difference between the results of the two tests was 0.84 mg/m?. The difference in direct measurements of the mass of 
the two samples should be no more than 0.1 mg. 

Discussion and Conclusion. The results obtained demonstrated the possibility of conducting a minimum number of 
measurements that, under reproducible conditions, are considered satisfactory and could be provided to experts for 
confirming the competence of the laboratory. As part of a special assessment of working conditions, employees of the 
object under study are recommended to use personal respiratory protection equipment that offers protection against 
highly dispersed dust particles. 


Keywords: convergence, intra-laboratory control, dust content, construction site, workplace, special assessment of 
working conditions (SAWC) 
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Введение. Измерения Ha запыленность рабочих мест проводятся B рамках специальной оценки условий 
труда (СОУТ) на строительной площадке с целью выявления и оценки вредных и опасных производственных 
факторов, а также разработки мероприятий по их устранению или минимизации [1]. Результаты специальной 
оценки условий труда являются основой для разработки и внедрения соответствующих мер по улучшению 
условий труда [2]. Это может быть изменение рабочих процессов, применение новых технологий, обновление 
оборудования или улучшение санитарных условий [3]. Важно, чтобы эти меры были реализованы с учетом 
мнения сотрудников и их потребностей [4]. 
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Проведение специальной оценки условий труда на строительной площадке имеет свои особенности, 
связанные с характером и условиями труда [5]. К таким особенностям можно отнести: 

— работу на открытом воздухе, что может привести к воздействию неблагоприятных погодных условий; 

— использование строительной техники и оборудования, которые могут создавать шум, вибрацию, пыль и 
другие загрязнения; 

— выполнение работ на высоте, что требует специального обучения и аттестации персонала; 

— работа с различными материалами и веществами, которые могут быть опасными для здоровья работников; 

— необходимость соблюдения требований охраны труда и техники безопасности, включая использование 
средств индивидуальной защиты. 

В результате выполнения COYT устанавливаются классы условий труда для каждой профессии, 
рекомендуются мероприятия по минимизации выявленных рисков и последствий влияния вредных 
факторов [6]. По данным СОУТ руководство организации разрабатывает мероприятия по улучшению условий и 
охраны труда, также определяются размеры дополнительных выплат для работников, занятых на производстве 
с неблагоприятными или рискованными условиями [7]. При инструментальном исследовании рабочих мест 
осуществляются замеры и анализ шума, степени запыленности, естественной и искусственной освещенности, 
температуры, влажности и пр. [8]. 

СОУТ имеют право проводить организации, обладающие соответствующей аккредитацией и опытом работы 
в данной области. Ими могут быть как специализированные организации, так и лаборатории, занимающиеся 
исследованием условий труда. 

В соответствии с Федеральным законом OT 28.12.2013 № 426 «О специальной оценке условий труда» 
организация, занимающаяся COYT, должна: 

— иметь не менее 5 штатных сотрудников с действующим аттестатом на право проведения СОУТ; 

— быть зарегистрированной в реестре организаций, проводящих СОУТ, который ведет Минтруда и 
соцзащиты РФ; 

— обладать системой менеджмента качества (СМК), соответствующей ГОСТ ISO/TEC 17025-2019 «Общие 
требования к компетентности испытательных и калибровочных лабораторий» и ГОСТР 54934-2012/ 
OHSAS 18001:2007 «Системы менеджмента безопасности труда и охраны здоровья. Требования»; 

— располагать аккредитованной испытательной лабораторией, которая выполняет исследования и измерения 
опасных и вредных факторов. 

Испытательная лаборатория (ИЛ) должна иметь действующий аттестат аккредитации, соответствующую 
область аккредитации, отражающую возможность исследования вредных и опасных факторов производства и 
условий труда, а также аттестованное испытательное оборудование (ИО) и поверенные средства 
измерений (СИ) [9]. Согласно постановлению Правительства РФ от 26 ноября 2021 года № 2050 «Об 
утверждении правил осуществления аккредитации в национальной системе аккредитации...» аккредитованные 
лаборатории в рамках подтверждения компетентности должны проводить сравнительные испытания для 
определения сходимости результатов измерений. Сходимость результатов измерений при проведении 
специальной оценки условий труда является одним из ключевых показателей, определяющих качество 
проводимых измерений. 

Проведенные ранее исследования в данной области затрагивали в основном межлабораторные 
сравнительные испытания (МСИ) и не рассматривали внутрилабораторные сравнительные (сличительные) 
испытания. Так, в работе Козлова А.В. приведены результаты МСИ геосинтетического материала [10]. Автор 
указывает на целесообразность нормирования требований показателей прецизионности испытаний. В 
исследовании Лекомцевой М.М. и Шендалевой Е.В. проведен анализ результатов МСИ нефтепродуктов, по 
результатам которого рекомендуется преобразовывать данные в линейную зависимость для получения 
неизменного значения стандартного отклонения [11]. В работе Ясина Дургута отмечается, что 
межлабораторные сравнительные испытания, проводимые в соответствии с ГОСТ ІЅОЛЕС 17025-2019 «Общие 
требования к компетентности испытательных и калибровочных лабораторий», выполняются более эффективно, 
по сравнению со сравнительными измерениями, организованными в соответствии с ГОСТ ІЅОЛЕС 17043 
«Оценка соответствия. Основные требования к проведению проверки квалификации» [12]. 

Между тем, согласно ГОСТ ІЅОЛЕС 17025-2019 «Общие требования к компетентности испытательных и 
калибровочных лабораторий» весьма важны еще и результаты внутрилабораторных сравнительных испытаний, 
так как они используются для проведения корректирующих мероприятий, улучшения качества работы 
лаборатории и повышения точности результатов [13]. Поэтому целью данной работы является оценка качества 
проведения выездных внутрилабораторных сличительных испытаний. Такая оценка представлена на примере 
определения запыленности рабочего места каменщика на строительной площадке. Лаборатории, 
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осуществляющие свою деятельность в рамках COYT, вынуждены проводить внутрилабораторные 
сравнительные испытания на выездной основе, и условия их значительно отличаются от лабораторных. 
Сотрудники лаборатории ограничены временем допуска на объект, поэтому им важно выполнить исследования 
не только качественно, но и с минимально количеством испытаний. В связи с этим задача авторов 
статьи — определить это минимально необходимое и достаточное количество таких испытаний. 

Методы исследований. Общие санитарно-гигиенические требования к воздуху рабочей зоны, а также 
методы контроля показателей микроклимата и содержания вредных веществ в воздухе регламентирует 
ГОСТ 12.1.005-88 «Межгосударственный стандарт. Система стандартов безопасности труда. Общие 
санитарно-гигиенические требования к воздуху рабочей зоны». Для определения запыленности использовался 
весовой метод!, который заключается в сборе пыли на фильтры с последующим взвешиванием и определением 
концентрации. 

Массовую концентрацию K, пыли определяют по формуле: 


(m, — mg ):1000 


Ky = у = , (1) 
20 
где К, — концентрация пыли в воздухе, мг/м?; то — масса чистого фильтра, MI; Mn — масса фильтра с 
осевшими частицами пыли, мг; V5o = объем воздуха, приведенный к стандартным условиям, ДМ?. 
V,.293.P 
Vy = (2) 


(273+Т)-101,33' 


где V, — объем воздуха, прошедший через фильтр, дм?; P — атмосферное давление, кПа; Т — температура 
воздуха на рабочем месте, °С. 

Согласно рекомендациям государственной системы обеспечения единства измерений? для обработки 
результатов определения запыленности рабочего места можно использовать метод приведения, который 
предполагает наличие значений измеряемых аргументов, полученных в результате многократных измерений. 

Результат косвенного измерения находят по формуле: 


* L A y 
A= ; — (3) 
jAL 
где L — число частных значений измеряемой величины; A; — j-e значение измеряемой величины; j — номер 


последовательного измерения от 1 до Г. 


Среднее квадратическое отклонение 5(4) случайных погрешностей результата косвенного измерения 


рассчитывают по следующей формуле: 


в (4-4) 
YE L(L-1) 9 
При нормальном распределении отдельных значений измеряемой величины доверительные границы 
случайных погрешностей определяют по формуле: 
A=t,-S(A), (5) 
где £j — коэффициент Стьюдента, зависящий от доверительной вероятности и числа результатов наблюдений. 


Прецизионность оценивается по показателю повторяемости (относительной среднеквадратического 
отклонения повторяемости), О, И по показателю точности (границ относительной погрешности), б: 


gs (6) 
где x; — 1-Й результат испытаний, полученный в условиях повторяемости; X — среднеарифметическое 
значение из п результатов испытаний в условиях повторяемости (сходимости). 

А 0, 0 
8 =—100% < 25%. (7) 
X 


1 Методика измерений массовой концентрации пыли гравиметрическим методом для целей специальной оценки условий труда. МИ 
АПФД-18.01.2018. URL: https://normativ.kontur.ru/document?moduleld- 1 &documentId-390372 (дата обращения: 25.11.2023). 

2 Рекомендация. Государственная система обеспечения единства измерений. Измерения косвенные. Определение результатов измерений и 
оценивание их погрешностей. МИ 2083-90. Москва: Комитет стандартизации и метрологии СССР; 1991. 
URL:https://znaytovar.ru/gost/2/MI 208390 GSI Izmereniya kosve.html (дата обращения: 25.11.2023). 
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В процентном отношении сходимость определяют в соответствии с выражением: 


Oo, 
r=—_ 100%. (8) 
X 
Для оценки сходимости результатов двух измерений необходимо использовать формулу: 
_ -æl ооо 
p, = 2 100 %, (9) 
X 


где X1, X2 — результаты двух параллельных измерений. 
Результаты измерений, полученные B условиях сходимости, признаются удовлетворительными, если 
полученная СХОДИМОСТЬ r будет меньше или равна нормативной CXO/IHMOCTH Гнорм: 
r € ку 1596. (10) 


норм 
За результат принимается выражение (8): 
х= Х+А. (11) 

Минимальное количество наблюдений зависит от коэффициента вариации, т. е. отношения 
среднеквадратического отклонения к среднему значению, и выбирается предварительно из опытных данных. 
Ошибку в определении среднего значения (статистическая погрешность измерений) необходимо искать в 
предположении, что закон распределения генеральной совокупности измерений — нормальный. Зная ошибку и 
доверительную вероятность, можно однозначно установить доверительный интервал для среднего значения и 
оценить максимальную статистическую погрешность измерений как отношение половины разности границ 
ошибки среднего значения к среднему значению из опытных данных [14]. 

Результаты исследований. Проведены измерения на определение запыленности строительной площадки 
при строительстве многоквартирного жилого дома. Отбор проб производился на рабочем месте каменщика в 
зоне кладочной стены фильтрами типа АФА (аналитические аэрозольные фильтры) с помощью двух 
аспираторов ПУ-4Э, обеспечивающих отбор проб с заданным объемным расходом через поглотитель по 
четырем параллельным каналам. Для установления статистической погрешности измерений инструментами и 
сотрудниками лаборатории при заданной доверительной вероятности Р = 0,95 выполнено восемь параллельных 
измерений, результаты которых приведены в таблице 1. 


Таблица 1 
Результаты восьми параллельных испытаний 
c a x 
E : = “> © Б ы 
= = Er m ш ES aes a o o р 
= = e i= i= 
Е £ Е E E B на шд 3% Е © 
S ij E Е | ва ЕЕЕ ЕН Ё | 
D Е, а. я Е z = 5 Ба ан чо S. E 
S E E E Я 8, fe йш ms Я = | > 
о SOS d rz Wo [5] 
z EB E e [5] > eo o я о = mo © © A 
e я я o o ral с 5 а Я о = я T 
ү а | В ЕЕ 391 |р 
2 З 5 Е | ТШЕ ДЕ 22/2 |$ 
© © x = = © © T [on z B E E E 
Es Es g 9 D om 7858 o © 5 Е 
> > À = К E = = а, © я 
Е ЕЯ 
< е) ©) 
п=8 
Nol 30965,6 30966,5 2 102,2 60 21,1 7,46 
No2 31346,8 31347,8 2 102,2 60 21,1 8,29 
№3 33980,5 33981,6 2 102,2 60 21,1 9,12 
№4 35572,3 35573,6 2 102,2 60 21,1 10,78 | 9,33 1,39 10,29 | 14,85 
№5 34926,9 34928,3 2 102,2 60 21,1 11,61 
№6 30477,7 30479,0 2 102,2 60 21,1 10,78 
№7 32678,4 32679,5 2 102,2 60 21,1 9,12 
№8 32678,4 32679,6 2 102,2 60 21,1 9,95 


Чтобы проверить, подчиняются ли данные нормальному закону распределения, построена гистограмма 
частот как один из способов визуального представления распределения данных (рис. 1). Гистограмма имеет 
форму колокола и близка к кривой нормального распределения (рис. 2). 
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Рис. 1. Гистограмма частот при и = 8 


Рассчитаны и на графике функции плотности нормального распределения приведены среднее значение и 
доверительный интервал проведенных измерений (рис. 2). 
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Puc. 2. Графики нормального распределения при п = 8: 
а — функция плотности нормального распределения; 6 — интегральная функция распределения 


За результат принимается среднее значение с указанием доверительного интервала: 
Kpg = 933 + 0,96 мг/м*. 
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Для определения минимального количества измерений найдены коэффициент вариации, средняя ошибка 
выборки и предельная ошибка выборки (рис. 3). Максимальная статистическая погрешность составила 0,26 мг/м?. 


hesccccccccccsccccccccccccccceccccccccsseccccccccccccccccos 


© R = 4.15 мг/м? 
10,0 Среднее значение e 


Кв=0.15 0,720.09 м/м — цп=+1,16мг/м? 


Запыленность, Mr/M2 
e 
© 


0 1 2 3 4 5 6 7 8 9 
Количество измерений, шт. 


Рис. 3. Размах результатов измерений 


Минимальное количество измерений определяется по формуле [11]: 


2.2 

t 
Пит = x , (12) 
Пит = 1,20 = 2. 


Согласно полученному значению минимального количества испытаний проведены два параллельных 
испытания по определению концентрации пыли на рабочем месте каменщика (таблица 2). 


Таблица 2 
Результаты двух параллельных испытаний 
с E 3 + = 
У Р 2 «| S ЕЕ | gs | „ Ж = 
ES d 2. ЗЕ 8 os БЕ 5 Е = 5] $e) ах 
5 E с? Б & E а EIRE SE. 89,5 24] 8a 
a. A чох >, о o он d 4 ш S| ш XS т р © o Б © 
El BS) :| FL SE $8| Se lh gg Sage) FEE) ЕВ ES 
© 8. я Ё © S e| 32| БЕ ЕРЕ gej RR SE, КЕЕ 
2 8 Ө ч | #4 > СӨ | = 8 HB ЕЕ жле ss 
* 8 2 єс жЕ | ЕВ ЗН | BE В| ОБ ов 
> 9 m Balm Ө O 6$ © E e 
e E x x 
n=2 
Nel | 34563,8 | 34565,0 2 102,2 60 22,0 9,98 
9,57 0,59 8,50 8,68 
Ne2 | 30654,3 | 30655,4 2 102,2 60 22,0 9,15 


Доверительный интервал результатов двух определений концентрации пыли при вероятности 0,95 
составил 0,81 мг/м? (рис. 4). 
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Рис. 4. Функция плотности нормального распределения при п = 2 
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Согласно условиям сходимости по ГОСТ Р ИСО 5725-1-2002 «Точность (правильность и прецизионность) 
методов и результатов измерений. Часть 1. Основные положения и определения» и ГОСТР 51672-2000 
«Метрологическое обеспечение испытаний продукции для целей подтверждения соответствия. Основные 
положения» результаты параллельных измерений можно признать удовлетворительными, так как выполняются 
условия (7) и (10) (рис. 5): 


r, = 8.68 % < поры = 15 %, 5, = 8,50 % < би = 25 96. 

30,0 30,0 
250 |---------------------------- 25,0 
20,0 20,0 


Сходимость результатов, % 
Относительная погрешность, % 


2 3 4 5 6 7 8 


Количество измерений, шт. 


Сходимость 

= = - Предельный показатель сходимости 

Относительная погрешность 

= - ~ Максимальная граница относительной погрешности 


Рис. 5. Зависимость сходимости результатов измерений и относительной погрешности 
от количества испытаний 


По результатам обработки принимается: 
Кир =9,57+0,81 мг/м?. 


Определение концентрации пыли весовым методом имеет достаточно высокую составляющую 
погрешностей используемых средств измерений и проводимых операций. Максимально допустимая разница 
результатов двух измерений составляет 0,84 мг/м?. Разница прямых измерений массы двух образцов должна 
быть не более 0,1 мг. При разнице масс в 0,2 мг показатель сходимости результатов двух параллельных 
измерений становится более 18 % и не удовлетворяет условиям прецизионности. При проведении испытаний на 
определение запыленности рассмотренной методикой возможно ограничиться двумя параллельными 
измерениями, при этом погрешность измерений будет равняться инструментальной погрешности. 

Строительная пыль относится к малоопасным (ТУ класс), разовая предельно допустимая 
концентрация (ПДК) составляет 6 мг/м?, суточная — 10 мг/м?. По результатам испытаний выявлено 
превышение разовой ПДК в 1,6 раза. Для снижения концентрации пыли на рабочем месте каменщика 
рекомендуются комплексные меры по уменьшению пылеобразования на строительной площадке, 
т.к. основными источниками пыли являются подъездные дороги, места выгрузки и погрузки строительных 
материалов, а также технологические процессы, связанные с обработкой и резкой строительных 
материалов [15]. Необходимо проводить мероприятия, снижающие концентрацию пыли на строительной 
площадке и за ее пределами [16]. 

Гигиенический норматив ГН 2.1.6.3492-17 «Предельно допустимые концентрации загрязняющих веществ в 
атмосферном воздухе городских и сельских поселений» устанавливает разовую ПДК строительной пыли в 
городском воздухе не более 0,5 мг/м?. 

Обсуждение и заключение. Определение концентрации пыли весовым методом имеет достаточно высокую 
составляющую погрешностей используемых средств измерений и проводимых операций. Выявлено, что 
максимально допустимая разница результатов двух проведенных испытаний составила 0,84 мг/м?. Разница 
прямых измерений массы двух образцов должна быть не более 0,1 мг. 

Обработка результатов испытаний показала возможность проведения двух измерений по определению 
запыленности рабочего места каменщика на строительной площадке. Сходимость результатов испытаний в 
рамках подтверждения компетентности лаборатории является важным показателем качества проведенной 
оценки и свидетельствует о достоверности полученных данных, поэтому в рамках внутрилабораторных 
сличительных испытаний рекомендуется использовать не менее шести измерений [17]. Чем выше сходимость 
результатов, тем меньше вероятность ошибок и неточностей в оценке условий труда, что, в свою очередь, 
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обеспечивает более объективную информацию о состоянии условий труда на рабочих местах и позволяет 
принять обоснованные меры по их улучшению [18]. 

Работникам исследованного объекта рекомендовано использовать средства индивидуальной защиты органов 
дыхания, защищающие от высокодисперсной пыли. Для снижения концентрации строительной пыли на 
строительной площадке, а также за её пределами следует проводить орошение строительной площадки водой 
для осаждения высокодисперсной пыли, поливать подъездные дороги, промывать колеса автотранспортных 
средств как на въезде на строительную площадку, так и на выезде [19]. 
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Аннотация 

Введение. Вредные компоненты рудной пыли, образующейся при разгрузке продуктов подготовки железорудного 
концентрата (ППЖК) на стадии смешивания, наносят ущерб как работникам, так и оборудованию. Для решения 
этой проблемы применяется орошение жидкостным аэрозолем с использованием форсунок больших диаметров 
(>20 мкм). Однако данный метод неэффективен в улавливании мелкодисперсных частиц пыли, поэтому 
повышение эффективности метода осаждения орошением пыли ППЖК становится актуальной задачей. Целью 
данного исследования является изучение воздействия технологии «Сухой туман», генерирующей капли жидкости 
размером до 20 мкм, на этапе разгрузки ППЖК горно-металлургического предприятия при осаждении взвешенной 
мелкодисперсной пыли. Основной задачей данного исследования являлась оценка эффективности и возможных 
преимуществ применения технологии «Сухой туман» для орошения пыли с последующим осаждением, поскольку 
к пыли ППЖК описанная выше технология ранее не применялась. 

Материалы и методы. Эксперимент по осаждению пыли ППЖК проводился в специально созданном 
лабораторном стенде. Посредством физического моделирования были получены параметры процесса 
осаждения. Далее полученные результаты подвергались анализу с точки зрения получения зависимости 
осаждения пыли с течением времени с учетом влияния технологии «Сухой туман». Для физического 
моделирования была разработана программа эксперимента. Согласно данной программе, пыль равномерно 
загружалась внутрь лабораторного стенда (сверху), распределялась в воздушном потоке по всему объему 
стенда крыльчаткой, а прибор, расположенный в нижней части, фиксировал изменение концентрации во 
времени. Далее были проведены эксперименты по осаждению пыли с применением жидкостного орошения. 
Совместно с загрузкой пыли в объем лабораторного стенда посредством форсунок, генерирующих капли 
размером 10 и 15 мкм, подавалась жидкость — отфильтрованная вода. Эффективность технологии «Сухой 
туман» при осаждении пыли ППЖК определялась визуально, и далее — на основании сопоставления графиков. 
Изучалась динамика изменения усредненных концентраций пыли от времени как при осаждении без орошения, 
так и с применением технологии «Сухой туман». В процессе эксперимента фиксировались характеристики 
микроклимата внутри лабораторного стенда (влажность, температура и скорость движения воздуха) и 
параметры двух форсунок — их рабочее давление и время распыления подаваемой жидкости. 

Результаты исследования. Сравнение результатов эксперимента показало уменьшение времени осаждения на 
40 % и 75 % при использовании форсунок на 10 мкм и 15 мкм соответственно. 

Обсуждение и заключение. По результатам эксперимента подтверждена эффективность технологии «Сухой 
туман» для осаждения пыли ППЖК при разгрузке на стадии смешивания. Полученные базисные результаты 
позволят в дальнейшем оценить эффективность осаждения пыли с применением дополнительно режима 
пульсирующей вентиляции. В таком сочетании ожидается повышение эффективности осаждения еще 
на 20-25 % относительно результатов, представленных в данной статье. Полученные результаты дают 
возможность обосновать рациональные параметры и применить на производстве вышеописанный способ для 
повышения эффективности осаждения пыли. Помимо этого, они создают основу для разработки методики 
ускорения осаждения пыли ППЖК с применением метода пульсирующей вентиляции. 
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Abstract 

Introduction. Harmful components of ore dust, formed during the unloading of products in the preparation of iron ore 
concentrate (PPIOC) at the mixing stage, cause damage to both workers and equipment. To address this issue, liquid 
aerosol spraying using nozzles with large diameters (>20 um) is used. However, this method proves ineffective in 
capturing fine-dust particles. Therefore, enhancing the efficiency of the dust deposition method through PPIOC dust 
spraying becomes a pressing challenge. The aim of this study is to investigate the impact of the Dry Fog technology, 
generating liquid droplets up to 20 um in size, during the unloading stage of PPIOC at a mining and metallurgical 
enterprise in the precipitation of suspended fine-dust particles. The primary goal of this research was to assess the 
effectiveness and potential advantages of applying the Dry Fog technology for dust spraying with subsequent 
precipitation, as this technology has not been previously applied to PPIOC dust. 

Materials and Methods. The experiment on the PPIOC dust deposition was conducted in a specially designed 
laboratory setup. Through physical modeling in the laboratory setup, parameters of the precipitation process were 
obtained. Subsequently, the results were analyzed to understand the dependence of dust precipitation over time, taking 
into account the influence of the Dry Fog technology. An experiment program was developed for physical modeling. 
According to the devised program, dust was uniformly loaded into the interior of the laboratory setup (from the top), 
distributed in the air stream throughout the volume of the setup by a fan, and an instrument located at the bottom 
recorded changes in concentration over time. Experiments on dust precipitation were then conducted using liquid 
spraying (filtered water as the liquid) introduced into the setup through nozzles generating droplets with sizes of 10 and 
15 um, concurrently with the loading of dust into the laboratory setup. The effectiveness of the Dry Fog technology in 
the deposition of PPIOC dust was determined visually and further analyzed based on a comparison of graphs. The 
dynamics of changes in the average dust concentrations depending on time was studied both during precipitation 
without spraying and using the Dry Fog technology. During the experiment, the characteristics of the microclimate 
inside the laboratory setup (humidity, temperature and air velocity) and the parameters of two nozzles — their operating 
pressure and the supplied liquid spraying time — were recorded. 

Results. The comparison of the results showed a reduction in the dust precipitation time by 40 % and 75 % when using 
nozzles with sizes of 10 um and 15 jum, respectively. 

Discussion and Conclusion. The experiment results confirm the effectiveness of the Dry Fog technology for PPIOC 
dust precipitation during unloading at the mixing stage. Fundamental findings have been obtained, providing a basis for 
further assessment of the efficiency of dust precipitation with the additional application of pulsating ventilation. In such 
a combination, an additional 20—25 % increase in precipitation efficiency is anticipated compared to the results 
presented in this article. The obtained results will support the justification of rational parameters and the implementation 
of the described method in production to enhance dust precipitation efficiency. Additionally, they will aid in developing 
a methodology to accelerate the PPIOC dust precipitation using the pulsating ventilation method. 


Keywords: ore dust, dust from iron ore concentrate preparation products, iron ore concentrate, liquid spraying, mass 
transfer, dust deposition experiment, liquid aerosol, dust aerosol, pulsating ventilation method 
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Введение. В целях снижения концентрации атмосферной пыли Ha участках производства в 
горно-металлургическом комплексе успешно применяются разнообразные инновационные методы и передовые 
технологии обеспыливания. Используют фильтрационные системы, циклоны и вентиляционные установки, 
специально разработанные для эффективного улавливания и удаления мельчайших частиц пыли [1—3]. Помимо 
этого, инженеры применяют методы жидкостного орошения, аэрофлотации и другие передовые технологии, 
которые совместно обеспечивают комплексный подход к решению проблемы запыленности. Эффективность 
предпринимаемых мер по обеспыливанию тесно связана с уникальными особенностями производственных 
процессов, характеристиками применяемого оборудования, конструктивными особенностями помещений и 
другими важнейшими факторами. Систематическая забота о техническом состоянии оборудования не просто 
является частью технологического процесса, а выступает важным элементом стратегии по снижению 
воздействия производственной деятельности на окружающую среду, способствует оптимизации 
производственной эффективности. Комплексное внедрение вышеперечисленных мероприятий — совмещение 
различных методов и технологий — позволяют достичь максимальной эффективности в уменьшении уровня 
пылевой нагрузки в пределах рабочей зоны. 

На этапе разгрузки и подготовки железорудного концентрата, где формируется значительное количество 
мельчайшей дисперсной пыли, продолжительное время используется система жидкостного орошения [4—6]. 
Для создания мельчайшего жидкостного аэрозоля применяются разнообразные виды форсунок, образующих 
капли диаметром от 30 до 150 мкм. Помимо этого, используются форсунки и с более крупным диаметром. 
Время оседания t мельчайшей дисперсной пыли (4 = 1-10 мкм) при использовании таких форсунок может быть 
достаточно продолжительным и достигать 25 400 с (около 7 часов) [6]. На других этапах производства для 
эффективного пылеосаждения применяется передовая технология — «Сухой туман». В этой технологии 
используются форсунки меньшего диаметра, распыляющие жидкость с дисперсностью капель в диапазоне 
от 1 до 20 мкм [7]. В данной технологии каждая капля жидкости служит эффективным инструментом для 
захвата и обеспечения оседания самых мельчайших частиц, создавая уникальное сочетание технологии и 
инженерного искусства в борьбе с проблемой атмосферной запыленности. 

Исследований, связанных с применением вышеописанной технологии для осаждения мелкодисперсной 
пыли при подготовке продуктов железорудного концентрата, ранее не проводилось. Поэтому целью данной 
работы был анализ воздействия технологии «Сухой туман», при которой генерируются капли жидкости 
размером до 20 мкм, на этапе разгрузки IIIDKK горно-металлургического предприятия при осаждении 
взвешенной мелкодисперсной пыли. Задачей стала оценка эффективности и возможных преимуществ 
применения технологии «Сухой туман» для орошения пыли с последующим осаждением. 

Материалы и методы. Для оценки эффективности инновационной технологии «Сухой туман» провели серию 
экспериментов пылеосаждения мельчайших частиц продуктов железорудного концентрата с использованием 
форсунок, имеющих диаметры 10 и 15 мкм. Для выполнения эксперимента был разработан лабораторный 
стенд (рис. 1). Данный стенд использовался для анализа осаждения угольной пыли в [8]. Он представляет собой 
кубическую емкость объемом 1 м?, изготовленную из органического стекла в алюминиевой раме. 

В роли средства для мониторинга концентрации частиц в атмосфере использовался высокоточный измеритель 
массовой концентрации аэрозольных частиц — АЭРОКОН-П. Данный измерительный прибор был разработан с 
учетом требований к определению массовой концентрации пыли с разнообразным происхождением и химическим 
составом. Измеритель массовой концентрации обладает уникальной способностью регистрировать дисперсность 
исследуемых частиц, включая те, которые имеют диаметр до 10 мкм. 


Техносферная безопасность 


https:// bps-journal.ru 


Курносов И.Ю. и др. Оценка влияния системы орошения на этапе подготовки продуктов железорудного концентрата 


АГП, Ш b I, II, Ш 


e B I, II, III 


Puc. 1. Лабораторный стенд для осаждения пыли co схемой размещения оборудования и точками измерений: 
1 — генератор воздушного потока; 2 — высокоточный измеритель массовой концентрации аэрозольных 
частиц — «АЭРОКОН-П»; 3 — генератор тумана; 4 — метеометр; 5 — датчик влажности; 6 — датчик температуры; 
7 — форсунка; 8 — засыпная воронка; 9 — воздушный пульсатор 


Порядок проведения измерений обозначены на рис.1 буквами А, Б, В (расположение датчика 
Аэрокон-П) (2). Значения І, П и Ш (находящиеся рядом с буквами) соответствуют точкам проведения 
измерений. Согласно вышеописанной схеме, эксперимент проводился по всей площади дна лабораторного 
бокса, а именно в 9-и его точках (точки измерения I, II, Ш следует понимать, как располагаемые в перспективе, 
т.е. друг за другом). 

Эксперименты по осаждению пыли производились при контроле параметров микроклимата. Данные 
параметры на момент проведения эксперимента по осаждению контролировались при помощи метеометра 
TESTO 435 (4). При проведении экспериментов по осаждению пыли метеометр использовался для контроля 
параметров пылевоздушной среды в лабораторном боксе — температуры и влажности воздуха. Это позволило 
установить начальные условия эксперимента и провести его в контролируемых условиях. Описанный выше 
измерительный прибор был задействован для мониторинга и поддержания влажности воздуха внутри стенда в 
пределах 25-30 % и температуры — в диапазоне от 22 °C до 25 °С. Оценка скорости движения воздушного 
потока, исходящего из генератора, проводилась с использованием цифрового крыльчатого анемометра, что 
обеспечивало стабильную скорость движения воздуха на уровне 4 м/с. Для моделирования системы орошения в 
лабораторном стенде использовалась установка Е218, предназначенная для туманообразования посредством 
распыления жидкости из различных форсунок (рис. 1). Рабочее давление данной установки составляет 5,4 МПа, 
а максимальное — 12,41 МПа. Для эксперимента в данной установке использовались форсунки диаметром 
15 мкм и 10 мкм (рис. 1). 

Размеры частиц пыли, используемой при проведении эксперимента по осаждению, согласно 
гранулометрического анализа, находились в диапазоне OT 1 мкм до 40 мкм (рис. 2). Поскольку прибор для 
регистрации концентрации частиц в воздухе измеряет частицы диаметром до 10 мкм, требовалось определить 
процентное содержание этих частиц в пробе пыли. Согласно данному анализу, требуемый размер частиц пыли, 
используемый для проведения эксперимента по осаждению, составляет примерно 10 % от 1 грамма навески, 
используемой при анализе пыли [9]. Для корректного определения прибором изменения концентрации навеска 
исследуемой пыли была увеличена до 5 грамм с целью повышения концентрации мелкодисперсной пыли с 
диаметром от | мкм до10 мкм. 

Эксперимент по осаждению пыли был реализован в лабораторном стенде, описанном ранее, с соблюдением 
всех параметров микроклимата, регулируемых соответствующими измерительными устройствами. Заданное 
количество пыли массой т=5г вводилось в объем лабораторного бокса в течение 3-5 секунд с 
использованием засыпной воронки (рис. 1). Во время этого процесса генератор воздушного потока был 
активирован, обеспечивая подачу воздуха со скоростью У = 4 м/с. Измерительное устройство АЭРОКОН-П, 
используемое в роли измерителя массовой концентрации аэрозольных частиц, периодически регистрировало 
данные о концентрации пыли с интервалом в 5 секунд, автоматически передавая полученные значения на 
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монитор компьютера [10]. Этап эксперимента по осаждению пыли завершался в момент достижения значения 
концентрации пыли, зарегистрированного прибором, равным лп = 0,00 мг/м?. Достижение данного значения 
считалось индикатором завершения эксперимента, при котором весь пылевой материал в лабораторном боксе 


считался осевшим. 
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Рис. 2. Гранулометрический анализ пыли агломерационного производства 
горно-металлургического комплекса: 
С — процентное распределение количества частиц в пробе по размерам, %; 
dC — интегральное распределение процентного количества частиц в пробе по размерам, %; 
d — размер частиц в пробе 


Используя аналогичные параметры температуры, влажности и скорости движения воздуха при проведении 
эксперимента по пылеосаждению, исследования осаждения пылеводяного аэрозоля проводились при 
поочередном использовании 2-х форсунок (d = 10 и d = 15 мкм) и контроле подаваемого давления (р = 5,4 МПа) 
жидкости. Помимо этого, осуществлялся контроль времени распыления t— 2 мин (по истечению 2 мин 
влажность воздуха в лабораторном стенде становится максимальной и составляет 98,5 %). 

Результаты исследования. В описанных условиях для обеспечения достоверности экспериментальных 
данных была выполнена серия из 10 экспериментов по осаждению пыли. Результаты данных экспериментов 
были обработаны, сведены в графики и проанализированы. Графики изменения значений концентрации пыли с 
течением времени при проведении 10 экспериментов от времени осаждения приведены на рис. 3. 

В соответствии с графиком, представленном на рис. 3, среднее время осаждения пыли при отсутствии 
орошения и использовании пульсирующей вентиляции составляет в среднем 1 828 секунды (30,5 минуты). В 
течение этого времени пыль в условиях производства находится в воздухе рабочей зоны и наносит вред не 
только оборудованию, но и персоналу предприятия. 

Результаты 10 экспериментов по осаждению пылеводяного аэрозоля с использованием 2-х форсунок также 
были представлены в виде графиков динамики изменения средних значений концентрации от времени 
осаждения аэрозолей (рис. 4, 5). 
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Puc. 3. Графики изменения концентрации пыли or времени по результатам 10 экспериментов 
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Согласно графикам, представленным на рис. 4, среднее значение времени осаждения пылеводяного аэрозоля 
при использовании форсунок диаметром 10 мкм, составляет 1115 с (18,5 мин). Эффективность данного способа, 
относительно времени самоосаждения пыли, составляет порядка 40%. Однако, согласно рис.5, при 
использовании форсунки диаметром 15 мкм, среднее время осаждения пылеводяного аэрозоля составляет 475 с 
(= 8 мин.). Эффективность данного способа составляет 74 %. 
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Рис. 4. Графики изменения концентрации пылеводяного аэрозоля от времени (диметр форсунки 10 мкм) 


* 30 
[an 
2 
g 25 
© 
ra 
+ 
E 20 
A 
Б, 
5 15 
8 
Q 
o 
S 10 
я 
я 
Е 
[ss] 
E 5 
о 
Е 
5 
x 0 


— —-2 з —4 =j 6 =] — о —0 


Рис. 5. Графики изменения концентрации пылеводяного аэрозоля от времени (диметр форсунки 15 мкм) 


На рис. 6 представлены пересчитанные зависимости концентрации аэрозолей, полученные в результате 
экспериментов по орошению с использованием форсунок диаметром 10 мкм и 15 мкм. В ходе экспериментов 
были обнаружены разные концентрации пыли, которые были зарегистрированы прибором. Для удобства 
представления эти концентрации были пересчитаны и выражены в виде процентного соотношения, зависящего 
от времени осаждения. 

При сравнении определенных ранее значений эффективности осаждения пыли с применением форсунок 
на 10 и 15 мкм, последние дают значительный эффект. С применением форсунки на 15 мкм значение 
эффективности осаждения пыли примерно в 2 раза выше, чем с форсункой на 10 мкм. При использовании 
форсунки на 15 мкм частицы пыли проявляют низкую склонность к слипанию между собой и менее 
эффективно формируют капли жидкости больших размеров, что снижает эффект «конденсата» в процессе 
орошения. Мелкодисперсные капли лучше заполняют пылевоздушную среду и могут использоваться в тех 
местах технологического производства, где капли больших размеров могут нанести ущерб. 


Безопасность техногенных и природных систем. 2024;8(1):41—49. eISSN 2541—9129 


100 

90 
< 80 
E 70 
É 60 
S 
= 50 
& 
E 40 
о 
Е 30 
o 
20 

10 

0 

0 200 400 600 800 1000 
Bpema, c 


Осаждение пыли орошением (форсунка 15 мкм) Осаждение пыли орошением (форсунка 10 мкм) 


Рис. 6. График изменения концентраций пыли (%) от времени в условиях пылеосаждения орошением при применении 
форсунок с диметрами 15 мкм и 10 мкм 


Обсуждение и заключение. В результате проведенных экспериментов по осаждению пыли в лабораторном 
стенде было установлено время осаждения пыли при двух различных сценариях: при естественном осаждении 
пыли и при осаждении пылеводяного аэрозоля (путем воздействия жидкостного аэрозоля с использованием 
форсунок диаметром 10 и 15 мкм). 

В дальнейшем для разработки методики по осаждению пыли необходимо использовать следующие 
параметры, как рациональные: 

— влажность воздуха внутри лабораторного — ф = 25—30 %; 

— температура воздуха B лабораторном стенде — T= 22—25 °С; 

— скорость движения воздуха, создаваемая генератором воздушного потока, должна быть на уровне V=4 м/с. 

Также следует использовать следующие технические параметры системы орошения: 

— для рассеивания жидкости рекомендуется использовать форсунки с диаметрами 10 мкм и 15 мкм; 

— рабочее давление в системе подачи жидкости — р = 5,4 МПа; 

— время распыления аэрозоля жидкости — 2 мин. 

Полученные результаты показывают разные значения эффективности осаждения пыли с применением 
форсунок, генерирующих разные диаметры капель, применяемых в процессе орошения. Описанные в статье 
форсунки демонстрируют значительное улучшение процесса пылеосаждения орошением при использовании 
форсунок на 10 и 15 мкм, нежели при использовании в процессе орошения форсунок, генерирующих капли 
размером больше 15 мкм [7]. Кроме этого, уменьшение количества влаги в воздушной среде производственных 
цехов будет оказывать меньшее негативное воздействие на оборудование. 

Анализ полученных данных подтвердил более высокую эффективность процесса осаждения пыли при 
использовании технологии «Сухой туман» во время разгрузки ППЖК на этапе смешивания. Полученные 
результаты предоставляют основу для дальнейшей оценки эффективности процесса осаждения пыли уже при 
использовании метода пульсирующей вентиляции. Ориентируясь на работы других авторов, использующих 
данную технологию [7], прогнозируется дополнительное повышение значений эффективности осаждения на 
20-25 %, относительно применения технологии «Сухой туман». Полученные в данной работе результаты будут 
использованы при разработке в дальнейшем методики осаждения пыли продуктов железорудного концентрата 
с использованием метода пульсирующей вентиляции. 

Полученные значения показателей эффективности технологии «Сухой туман» дают основание 
рекомендовать его для внедрения на производстве для осаждения пыли на участках горно-металлургических 
производств с высокой пылевой нагрузкой. 
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Обеспечение безопасных условий труда при воздействии 
на работников поездных бригад виброакустических факторов 
Д.А. Coxonon®, Е.И. Головина (5a 


Воронежский государственный технический университет, г. Воронеж, Российская Федерация 
u00111@vgasu.vrn.ru EDN: WMKRZO 


Аннотация 

Введение. Повышенная интенсивность шума и вибрации на железнодорожном транспорте может стать не только 
актуальной проблемой в области охраны труда, но и серьезной экологической проблемой в ближайшем будущем 
из-за необходимости постоянного увеличения веса грузовых поездов, повышения скорости движения 
пассажирских поездов и возрастания пропускной способности железнодорожных линий. В последнее время 
немало сделано для минимизации негативного воздействия шума и вибрации на железнодорожном транспорте. 
Однако проблема превышения показателей виброшумовых факторов остается актуальной. Уменьшение вредного 
воздействия шума и вибрации способствует улучшению условий труда работников поездных бригад и 
повышению комфортных условий для пассажиров. Как правило, основное внимание при анализе воздействия 
виброакустических факторов уделяется внешнему их воздействию на селитебную территорию и 
меньше — воздействию на подвижной состав. К тому же такие исследования не являются комплексными, так как 
не дают полного представления о воздействии шума и вибрации непосредственно на состав в движении. Цель 
данного исследования — анализ воздействия виброакустических факторов на работников поездных бригад и 
пассажиров для разработки метода их комплексной защиты, основанного на применении виброгасящих 
материалов. 

Материалы и методы. В ходе работе были изучены нормативные документы, проведен комплексный анализ 
актуальной информации по данной теме, использованы методики по расчету вибрации и виброускорения. 
Значения уровней звукового давления были получены с помощью шумомера SPM-101. Объектом исследования 
был выбран вагон типа «плацкарт» во время движения. В рамках исследования было измерено звуковое давление 
на выбранных участках дороги и определен класс условий труда. 

Результаты исследований. Результат анализа воздействия виброакустических факторов на работников 
поездных бригад свидетельствует в необходимости усиления комплексных мер их защиты. Полученный расчет 
силы вибрации стал основой предложенного авторами метода минимизации вредного воздействия 
виброакустических факторов, в основе которого использование в конструкции вагона вибро- и 
звукоизолирующих материалов, подходящих по необходимым акустическим параметрам. В частности, 
предложено с целью шумоизоляции на пол вагона уложить плотную резину, а на потолок — материал из 
неорганического волокна. В работе также дана экономическая оценка эффективности мероприятий по 
минимизации вредного виброшумового воздействия. 

Обсуждение и заключение. Предложенные методы минимизации шумового и вибрационного воздействия могут 
способствовать снижению уровня звукового давления и показателей вибрации до нормативных значений, что 
сократит негативное воздействие виброакустических факторов на работников поездной бригады и пассажиров. 
Рассмотренный комплекс по минимизации вибрации и шума включает в себя совокупность методов, в которых 
нашли применение бюджетные материалы, созданные на основе вторично переработанного сырья. 


Ключевые слова: звуковое давление, вибрация, шумовое загрязнение, шумоизоляция 
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Abstract 

Introduction. The increased intensity of noise and vibration of railway traffic may become not only an urgent problem in 
the field of occupational safety, but also a serious environmental problem in the near future due to the constant need to 
increase the weight of freight trains, the speed of passenger trains and the capacity of railway lines. Recently, a lot has 
been done to minimize the negative effects of noise and vibration in railway transport. However, the problem of exceeding 
the indicators of vibration noise factors remains relevant. Reducing the harmful effects helps to improve the working 
conditions of train crew employees and improve the comfortable conditions of transport passengers. As a rule, the main 
focus of the analysis is on the external impact of vibroacoustic factors on residential areas and less attention is paid to the 
impact on rolling stock. However, these studies do not provide a complete picture of how noise and vibration actually 
affect the train in motion. The aim of the study was to obtain the result of an analysis of the impact of vibroacoustic factors 
on train crew workers and to propose a method of comprehensive protection based on the use of vibration damping 
materials. 

Materials and Methods. \n the course of the work, regulatory documents were studied, a comprehensive analysis of 
relevant information on this topic was carried out, and methods for calculating vibration and vibration acceleration were 
used. The values of the sound pressure levels were obtained using a SPM-101 sound level meter. The object of the study 
was a carriage of the "reserved seat" type in the process of movement. As part of the research plan, sound pressure was 
measured at selected sites and a class of working conditions was determined. 

Results. The result of the analysis of the impact of vibroacoustic factors on train crew workers indicated the need to 
strengthen comprehensive measures to protect them. The obtained calculation of the vibration force became the basis of 
the method proposed by the authors to minimize the harmful effects of vibroacoustic factors, which was based on the use 
of vibration and sound insulation materials suitable for the necessary acoustic parameters in the construction of the car. 
In particular, it was proposed to cover the floor of the car with dense rubber and the ceiling of the car with an inorganic 
fiber material for sound insulation. The work also provides an economic assessment of the effectiveness of measures to 
minimize harmful vibration noise effects. 

Discussion and Conclusion. The proposed methods for minimizing noise and vibration impacts can help reduce the level 
of sound pressure and vibration indicators to standard values, which will significantly reduce the negative impact of 
vibroacoustic factors on train crew workers and railway passengers. The considered complex for minimizing vibration 
and noise includes a set of methods in which affordable materials based on recycled substances have found application. 
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Введение. Железнодорожный транспорт всегда вызывал и вызывает шум и вибрацию, которые B большинстве 
случаев не могут быть полностью устранены и, скорее всего, будут негативно влиять на окружающую среду 
невзирая на применение различных методов смягчения их воздействия. Проблема уменьшения влияния шума и 
вибрации железнодорожного транспорта на людей и воздействия на инженерные конструкции очень сложна из- 
за большого количества параметров и факторов распределения шумовой и вибрационной энергии. Негативное 
влияние виброакустических факторов растет с каждым годом в связи с увеличением количества и веса грузовых 
вагонов, большого их износа, несвоевременным техническим обслуживанием пассажирских вагонов, а также 
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нехваткой свободного пространства в городских районах, в результате чего новые деловые и жилые объекты 
размещаются в непосредственной близости от железной дороги. Проблема шума и вибрации от 
железнодорожного транспорта особенно остро стоит в крупных городах с большим железнодорожным 
сообщением главным образом из-за малого расстояния от железнодорожных путей до ближайших зданий. Этот 
факт ограничивает возможности для принятия мер по снижению уровня шума и вибрации. На сегодняшний день 
первостепенной задачей является защита от превышения звукового давления по системе «внешний источник 
шума — окружающая среда». Система подразумевает воздействие шума на близко расположенные объекты 
инфраструктуры, жилые зоны, природу. Защита от такого рода воздействий достигается установкой 
акустических экранов. С их помощью возможно значительно сократить звуковое давление. 

Малоизученной до сих пор является проблема комплексной защиты от виброакустического воздействия. 
Стоит также отметить и недостаток вибродемпфирующих устройств, которые были бы способны снижать 
совокупное воздействие шума и вибрации. По мнению авторов, научная новизна данного исследования 
заключается именно в том, что рассмотрен комплексный подход к решению проблемы минимизации воздействия 
виброакустических факторов и предложены простые шумопоглощающие материалы для массовой их установки. 

Вибрация от рельсового транспорта может оказывать негативные воздействия непосредственно или в 
сочетании с фоновым шумом. Значительное превышение нормативов уровня шума доставляет дискомфорт не 
только работникам и пассажирам, но и жителям прилегающих к железнодорожному полотну территорий, часто 
является значимой причиной ухудшения самочувствия и снижения работоспособности. Специфика шумового 
загрязнения на железнодорожном транспорте такова, что при долгосрочном воздействии в совокупности с 
вибрацией оно очень негативно влияет на работников ж/д транспорта [1]. Длительное и систематическое 
воздействие излишнего шумового давления может привести к развитию хронических заболеваний нервной 
системы, опорно-двигательного аппарата, быть причиной потери бдительности, ухудшения самочувствия, 
ощущения дискомфорта, чувства раздражения [2]. Постоянное воздействие шума даже на номинальных 
значениях (особенно ночью) может стать причиной нарушений сна и хронических неврозов [3]. Количество 
людей с нарушениями сна очень велико — большинство сталкиваются с воздействием постоянного шума в 50— 
60 дБ именно ночью, что при отсутствии других дневных источников шума является одним из основных 
факторов бессонницы [4]. 

Из всех вредных экологических факторов вибрация является одним из самых массовых. Под ее воздействием 
находится от 50 до 70 % населения [5]. Анализ вибраций, возникающих в результате движения подвижного 
состава, очень важен, так как железная дорога — основная транспортная система, охватывающая большую часть 
территории нашей страны и стран ближнего зарубежья. С каждым годом растут объемы и количество 
перевозимых по железной дороге грузов [6]. 

Рабочая вибрация рельсов возникает на стыке колеса и рельса (система «колесо — рельс») и распространяется 
не только через систему опор пути на почвенный покров и окружающие здания, но и непосредственно на 
подвижной состав. В некоторых случаях пассажиры могут непосредственно ощущать вибрацию, которую обычно 
называют вибрацией от взаимодействия подвижного состава и рельсового полотна. Колебания грунта, вызванные 
железнодорожным движением, в основном возникают в зоне контакта между колесом и рельсом [7]. Поэтому 
важно, чтобы действие этого механизма возникновения негативного фактора было сведено к минимуму. Основным 
источником возникновения вибрации на железной дороге является вибрация при качении, вызванная неидеальной 
рабочей поверхностью колеса и дефектами на рабочей поверхности рельса. Эти вибрации передаются колесу и 
путевым конструкциям, что приводит к возникновению превышения пороговых значений по основным октавным 
уровням. Как следствие, вибрация, передаваясь через колесные пары, тележку, вызывает вибрацию кузова вагона, 
что особо ощутимо на высокой скорости поезда. Вибрации высокой амплитуды могут привести к критическому 
состоянию узлов вагона, к урону конструкции рельсового подвижного состава. С целью сохранения оборудования 
вибрации такого рода следует сокращать не только в системе «колесо — рельс», но и непосредственно принимая 
меры по защите основных частей вагона, это, в свою очередь, минимизирует появление вторичной вибрации от 
корпуса вагона. Чтобы свести к минимуму механическое воздействие колеса и рельса, поверхность рельса должна 
быть как можно более гладкой. Вибрации, создаваемые конструкцией вагона и тележки, также должны быть 
сведены к минимуму, чтобы уменьшить количество вибраций, передаваемых на людей в вагоне [8]. Для создания 
системы с желаемыми акустическими и структурно-динамическими свойствами при низком воздействии на 
человека шума и вибрации необходимо учитывать конструкцию пути, колес, тележек, вагонов. Другим важным 
аспектом является обеспечение технического обслуживания этих деталей соответствующим образом 
(своевременное техническое обслуживание и плановый ремонт). 
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Актуальные методы снижения вибрации, как правило, делятся на две основные категории: категория 
пассивных мероприятий и категория активных мероприятий. В настоящее время чаще всего используются 
пассивные мероприятия. К ним относится применение традиционных виброгасителей, систем амортизации 
ударов разного рода (соударение вагонов) и изоляции основания [5]. Однако современные методы снижения 
вибрации должны быть направлены на пассивно-активное управление вибрацией, которое заключается в 
активном приложении силы, равной противоположным силам, создаваемым внешней вибрацией. Иногда для 
гашения вносимых вибраций предлагается использовать разработанные системы управления и гашения сил, 
которые создают вибрацию, а также конструкции, созданные на основе новых композитных материалов. В 
частности, можно использовать пьезоэлектрические элементы в качестве гашения излишней вибрации!. Суть 
этого метода заключается в том, что пьезоэлектрический элемент вырабатывает электрическую энергию при 
механическом растяжении и реагирует по форме деформации, если на него подается электрическая энергия. Еще 
одним методом минимизации шума и вибрации является метод применения дисковых тормозов вместо чугунных 
колодок. Альтернативой же более сложным в адаптации, установке и дорогим дисковым тормозам является 
замена чугунных блоков на агломерационные или композитные. Все названные выше методы направлены на 
основные узлы вагона. Авторами был проведен расчет виброизоляции в самом вагоне на пассажирских местах. 

Материалы и методы. В качестве материала исследования был рассмотрен вагон типа «плацкарт». Замеры 
проводились с помощью прибора для измерения звукового давления — шумомера SPM-101. Участки с наиболее 
негативным воздействием виброакустического фактора были выбраны исходя из конструктивных особенностей 
тележки вагона и строения внутренней структуры пространства вагона. Следовательно, план исследования 


включал в себя измерение звукового давления в вагоне и определение класса условий труда. 

Для анализа возможности снижения шума в источнике вагона типа «плацкарт» рассмотрена схема вагона и 
отмечены на ней участки наиболее негативного воздействия виброакустических факторов (рис. 1). На этих 
участках подразумевается наибольшее звуковое давление в вагоне и как следствие — критическое негативное 
воздействие виброакустических факторов при различных условиях (скорость движения, состояние 
железнодорожного полотна и др.). 
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Рис. 1. Значения звукового давления в вагоне по участкам (1, 2, 3) 


Следующим методом исследования минимизации вибрационного воздействия был расчет виброизоляции 
пассажирских мест в вагоне на основе методики?. Были применены виброзащитные вставки из резинового 
материала. Пол вагона у основания кресла вибрирует с частотой f= 50 Гц, с виброскоростью у = 0,4 м/с и массой 
места 22 кг. Примем массу пассажира в 70 кг. Модуль упругости: соп = 0,04 МПа, Е = 2.5 МПа, й = 0,1 м, 
коэффициент сопротивления & = 600 Н с/м. 

Определим массу, приходящуюся на вставки, когда пассажир будет сидеть. Масса пассажира, приходящаяся 
на сидение, составляет 70 %, значит: 

т. = 70.0,7 = 49 кг (1) 
тоб =49+ 22 = 71кг (2) 


Определим максимальный статический прогиб вставок: 


2= пзп 94. 5 216.49. (3) 
E 250 


д 


' Снижение уровня шума в кривых. Железные дороги мира. 2009;6:70-76. URL: https://zdmira.com/images/pdf/dm2009-06 70-76.pdf (дата 
обращения: 14.11.2023) 

? Пособие по акустической виброизоляции центробежных машин. Москва: Издательство литературы по строительству; 1973. 35 с. 
URL: https://meganorm.ru/Data2/1/4293801/4293801338.pdf (дата обращения: 14.11.2023) 
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Частотная характеристика собственных колебаний амортизированного пассажирского места будет: 


05 = 2.314. 2 _ =79¢!, 


NE 16-107 


о=2т 


Коэффициент передачи вибрации определяется по формуле: 


Найдем относительное демпфирование D и круговую частоту o: 


Получим: 


D= E _ 600 _ 0,05. 
209 -ть 2.79.71 
лы 

00 79 
1+(2-0,05-4)" 
iz = 0,071. 


П-у] +(2:0,05:4)? 


Затем определим вибрационную скорость на сидении: 


у = 0,4.0,071= 0,028 м/с. 


Частота колебания места составит: 


oO 79 
R= 


= — = = 12,57 Гц. 
2n 2-3,14 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


Результаты исследования. Согласно ГОСТу 33787-2019 для частоты 12,57 Гц виброскорость не должна 
превышать 0,0056 w/c?. Полученный в результате расчета показатель не соответствует данному требованию. Для 
его улучшения можно применить виброзащитный материал с другими характеристиками — таким может стать 
композитный лист из плотных резин (рис. 2). Благодаря более высокому коэффициенту сопротивления можно 


добиться минимизации вибрационного воздействия. 


ГОСТ 33787-2019. 


Рис. 2. Вариант расположения виброзащитного материала: 
1 — пассажирское сиденье; 2 — виброзащитная вставка 


Оборудование железнодорожного подвижного состава. Испытания на удар 


https://docs.cntd.ru/document/1200170805 (дата обращения: 14.12.2023). 
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В таблице 1 представлены результаты измерения уровней звукового давления на участках, представленных 
на рис. 1. Сравнив их с нормами уровня звука в вагоне по СП 2.5.3650-20, можно сделать вывод, что наблюдается 
превышение показателей“. 

Таблица 1 
Показатели шумового давления в вагоне по зонам 


Участок, № Уровень звука, дБ 
1 77,9 
2 74,1 
3 78,5 


Проведенные в данной работе расчеты позволили авторам сделать вывод о необходимости применения 
метода минимизации виброакустических факторов, основанного на использовании виброгасящих материалов. 
Его эффективность подтверждается представленным исследованием. Авторами предложен комплекс 
мероприятий, включающий в себя установку в обшивку потолка вагона в наиболее подверженных шуму зонах 
шумопоглощающего материала для минимизации шумового загрязнения. Материал представляет собой 
композит из неорганических волокон толщиной листа OT 2 до 4см. Также для снижения шумового и 
вибрационного воздействия предложено установить вставки в пол вагона по аналогии с созданием 
амортизирующей системы. Для расчета эффективности метода был проведен расчет воздействия вибрации на 
пассажирское место с учетом вибрации пола вагона. Материал безопасен при эксплуатации для человека и 
окружающей среды, соответствует критериям пожарной безопасности, которые установлены для пассажирского 
вагона согласно распоряжению ОАО «РЖД» от 5 ноября 2009 г. № 2255p°. 

Обсуждение и заключение. Оценка воздействия вибрации и шума, полученная в данном исследовании, 
показывает, что оно существенно влияет на условия труда работников поездных бригад и пассажиров. Работники, 
которые много времени проводят в дороге, рабочие, занимающиеся ремонтом путей в непосредственной 
близости от железнодорожного полотна, и пассажиры поездов подвергаются крайне вредному шумо- и 
вибрационному воздействию, которое оказывает разрушительное воздействие на их здоровье. Поэтому так важно 
минимизировать это вредное воздействие в источнике его возникновения. Несомненно, улучшение 
существующих и внедрение новых методик минимизации последствий шумового и вибрационного факторов на 
железной дороге положительно повлияет на здоровье работников и будет способствовать созданию 
благоприятной социальной обстановки. Стоит особо подчеркнуть, что мероприятия по минимизации негативного 
воздействия шума и вибрации от железнодорожных составов дают относительно небольшой экономический 
эффект, но имеют при этом большую ценность для обеспечения защиты здоровья как работников, так и 
пассажиров, создания для них более комфортных условий. 

В предложенном авторами комплексе мероприятий по минимизации вибрации и шума применяются 
материалы, созданные из вторично переработанного сырья. Для пола вагона используется плотная резина, для 
потолка — композит из неорганического волокна толщиной 2—4 см. Оснащением вагонов этими материалами 
можно достичь: 

— минимизации шумового давления в вагоне на 7-12 Дб; 

— уменьшения вибрационного воздействия. 
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Аннотация 

Введение. Сера присутствует в органической и минеральной части всех видов твердого топлива, причем от ее 
содержания зависит проявление опасных свойств шахтопластов. В связи с этим не теряют актуальность научные 
изыскания, связанные с определением количества серы в ископаемом топливе. Содержание общей серы и ее 
разновидностей изучалось в большинстве случаев с позиций эффективности процессов переработки топлива, качества 
и экологической безопасности угольной продукции. Меньше внимания уделяли проявлению опасных свойств 
шахтопластов, прогнозированию и разработке профилактических мероприятий по предотвращению эндогенных 
пожаров в угольных шахтах. Цель представленной работы — создание инженерного метода расчетного определения 
элементного содержания серы в минеральной массе углей для прогноза опасных свойств шахтопластов. 
Материалы и методы. По результатам корреляционного анализа установили, каким образом содержание серы 
в горючей части топлива зависит от общей серы. Рассмотрели данные по углям практически всех шахтопластов 
Донецкого и Львовско-Волынского бассейнов. Представили соотношения зольности пластовых и обогащенных 
проб, а также сведения о содержании общей серы и серы в горючей части. Это позволило расчетным путем 
оценить содержание серы в минеральных примесях. Исходные данные для разработки метода взяли из 
справочно-нормативных документов. 

Результаты исследования. Рассмотрены семь совокупностей, которые включали от 149 до 1827 шахтопластов. 
По каждой обобщили справочные и расчетные данные: 

— содержание серы — общей (S) и в горючей части топлива (5,); 

— условные точки для корректировки показателей во избежание получения отрицательных усредненных значений; 

— эмпирические уравнения, которые описывают расчетные и скорректированные зависимости по соотношению 
SP и S; 

— корреляцию этих показателей; 

— диапазоны изменения выхода золы. 

Все эти сведения представлены для обогащенных и не обогащенных углей, для выхода золы более и менее 10 %. 
В одном случае отмечены удовлетворительные результаты определения расчетного содержания серы в условно 
негорючей массе. Это совокупность из 149 шахтопластов. Для 13 из них разность показателей (5 — S) больше 
или равна нулю. Для 136 — меньше нуля. Для 81 содержание серы в условно негорючей части топлива $, меньше 
нуля. При этом отмечается достаточный уровень корреляции (0,79) и ограниченный диапазон изменения выхода 
золы (2,410). Соответствующее этой совокупности эмпирическое уравнение рекомендуется для прогноза 
опасных свойств шахтопластов с выходом золы менее 10 %. 

Обсуждение и заключение. Полученные результаты можно использовать для совершенствования нормативных 
документов по безопасному ведению горных работ с учетом содержания серы, как в органической массе, так и в 
минеральных примесях ископаемого угля. 


© Пронская H.B., Филатьева Э.Н., Филатьев M.B., Шашло H.B., 2024 


Пронская Н.В. и др. Onpedenenue содержания серы 6 минеральной массе для прогноза опасных свойств угольных шахтопластов 


Ключевые слова: сера в ископаемом угле, инженерный метод определения серы в минеральной массе угля, 
предотвращение эндогенных пожаров в угольных шахтах 


Благодарности. Авторы выражают благодарность редакционной коллегии журнала и рецензенту за 
профессиональный анализ и рекомендации по корректировке статьи. 


Для цитирования. Пронская Н.В., Филатьева Э.Н., Филатьев M.B., Шашло Н.В. Определение содержания серы 
в минеральной массе для прогноза опасных свойств угольных шахтопластов. Безопасность техногенных и 


природных систем. 2024;8(1):58—72. https://doi.org/10.23947/2541-9129-2024-8- 1-58-72 


Original article 


Determination of Sulfur Content in Mineral Mass for Prediction of Hazardous Properties 
of Coal Mine Seams 


Natalya V. Pronskaya ! ©, Elvira М. Filatieva ? ©, Mikhail V. Filatiev ? ©, Nina V. Shashlo ? (954 
! Donbass State Technical University, Alchevsk, Russian Federation 

? Lugansk Vladimir Dahl State University, Lugansk, Russian Federation 

3 Don State Technical University, Rostov-on-Don, Russian Federation 

ninelllsss@gmail.com 


Abstract 

Introduction. Sulfur is found in the organic and mineral parts of all types of solid fuels, and its concentration determines 
the extent of the hazardous characteristics of mine seams. Therefore, research on the determination of sulfur content in 
fossil fuels has not lost its relevance. The total sulfur content and its varieties has been studied in most cases in terms of 
the efficiency of fuel processing processes, the quality and environmental safety of coal products. Less attention was paid 
to the hazardous properties of mine seams, forecasting and developing preventive measures for endogenous fires in coal 
mines. This work aims to develop an engineering method to calculate the elemental sulfur content in the mineral mass of 
coal to predict the hazardous properties of mine seams. 

Materials and Methods. According to the results of the correlation analysis, we established how the sulfur content in the 
combustible part of the fuel depended on the total sulfur. We analyzed data on coals from almost all mine seams of the Donetsk 
and Lviv-Volyn coal basins. We presented the ash ratios for the formation and basin samples, as well as the information about 
the total sulfur and sulfur content in the combustible part. This allowed us to estimate the sulfur content in mineral impurities 
by calculation. The initial data for developing the method were taken from reference and regulatory documents. 

Results. Seven sets were considered, which included from 149 to 1827 mine seams. For each, reference and calculation 
data were summarized: 

— sulfur content — total (5/2) and in the combustible part of the fuel (5\); 

— conditional points for adjusting the indicators to avoid obtaining negative average values; 

— empirical equations that describe the calculated and corrected dependencies based on the ratio of S/ and Sr; 

— correlation of these indicators; 

— ranges of change in ash yield. 

All this information was presented for clean and raw coals, for ash yields of more and less than 10%. In one case, 
satisfactory results were noted for determining the calculated sulfur content in a conditionally non-combustible mass. 
This was a set of 149 mine layers. For 13 of them, the difference in indicators (S; — S,) was greater than or equal to zero. 
For 136 — less than zero. For 81, the sulfur content in the conditionally non-combustible part of the fuel 5, was less than 
zero. At the same time, there was a sufficient level of correlation (0.79) and a limited range of changes in ash yield 
(2.4—10). The empirical equation corresponding to this set was recommended for predicting the hazardous properties of 
mine seams with an ash yield of less than 10%. 

Discussion and Conclusion. The results of this study can be used to improve regulations for the safe conduct of mining 
activities, taking into account sulfur content in organic matter and in mineral impurities of fossil coal. 


Keywords: sulfur in fossil coal, engineering method for determining sulfur in the mineral mass of coal, prevention of 
endogenous fires in coal mines 
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Введение. Сера — элемент соединений, входящих в органическую и минеральную часть любого топлива. B 
минеральной массе углей сера встречается главным образом в виде сульфатов (в основном кальция и железа) и 
дисульфидов железа Fes (пирита и марказита), в органической массе — в виде органических соединений. 

Различают несколько основных форм серы в углях. Это части общей серы: 

— сульфатная Sso, — входит в состав неорганической массы угля в виде сульфатов металлов; 


— пиритная S, — входит в состав неорганической массы угля в виде дисульфидов металлов (пирита и марказита); 
— органическая S, — входит в состав органических соединений. 
Общую серу и ее разновидности, в отличие от серы золы, условно пересчитывают на элементную серу [1]: 

S, = 50, tSp +58, - (1) 


Содержание серы в золе принято рассчитывать в виде триоксида серы (SO3). Сера — нежелательная и вредная 
часть топлива. При сжигании угля она выделяется в виде SO», загрязняет окружающую среду, вызывает коррозию 
металлических поверхностей, уменьшает теплоту сгорания топлива. При коксовании сера переходит в кокс, 
ухудшая его свойства и качество металла. По указанным причинам содержание общей серы и ее разновидностей 
изучалось в большинстве случаев с позиций эффективности процессов переработки топлива, качества и 
экологической безопасности угольной продукции. В этом направлении получены хорошие результаты, что 
подтверждается разработкой и успешным использованием современной промышленной классификации углей по 
генетическим и технологическим параметрам. В значительной степени она базируется на данных справочников, 
в которых по итогам экспериментов обобщалась информация о качестве каменных углей и антрацитов [2], а 
также об их обогатимости [3]. 

Менее значимы успехи в прогнозировании и разработке профилактических мероприятий по предотвращению 
эндогенных пожаров при ведении работ в угольных шахтах. Так, в течение двадцати лет (с 1994 по 2015 год), в 
шахтах Украины произошло более двухсот эндогенных пожаров [4]. Согласно открытым данным, аналогичные 
аварии происходят во многих угледобывающих странах. Таким образом, не решена проблема, актуальная с 
XIX века. Вначале она была связана с самовозгоранием угля на складах и при транспортировке, затем стали 
происходить аварии в угольных шахтах, обусловленные эндогенными пожарами [5]. 

Содержание серы в угле не только определяет его потребительские качества, но и существенно влияет на 
проявление опасных свойств шахтопластов при ведении горных работ [6]. В частности, при увеличении 
содержания в угле общей серы на 1—6 % вероятность эндогенного пожара возрастает более чем в три раза [7]. 
Были попытки разобраться, как по минеральным примесям можно выявить опасные свойства угольных 
шахтопластов [8]. Однако и сейчас почти не известно, каким образом исследуемые риски зависят от вида и 
формы серы в ископаемом угле. В результате в нормативных документах l? вероятность эндогенных пожаров в 
шахтопластах определялась только содержанием общей серы?. 

Исследования [5] показали, что самовозгорание в значительной мере обусловлено природой не только угля, но 
и сопровождающих пород“. Крайне важно установить составную часть топлива, от которой зависит вероятность 
самовозгорания. Без этого лабораторные методы исследований не дают приемлемых ответов на поставленные 
вопросы. Предполагалось наряду с лабораторными исследованиями угольного вещества учитывать наличие в 
ископаемом топливе углистой породы из минеральных пропластков, а также из пород кровли и почвы. Для этого 
необходимо дополнительно рассматривать выход золы и содержание всех видов серы [5]. 

В руководстве по предупреждению и тушению эндогенных пожаров°, как и в других нормативных 
документах по безопасному ведению горных работ, не принимается во внимание наличие минеральных 
примесей. Они не учитываются при борьбе с пылью в угольных шахтах, прогнозах газовыделенияё и внезапных 
выбросов”? и в целом в нормативных документах, регламентирующих горные работы. 


! Руководство по предупреждению и тушению эндогенных пожаров на угольных шахтах Украины. КД 12.01.402-2000 от 18.12.2000. Донецк: 
Научно-исследовательский институт горноспасательного дела; 2000. 216 с. 

? Каталог углей СССР, склонных к самовозгоранию. Москва: Недра; 1982. 416 c. 

3 Руководство по предупреждению и тушению эндогенных пожаров на угольных шахтах Украины. КД 12.01.402-2000. Донецк: НИИГД; 2000. 216 с. 
^ ГОСТ 25543-2013. Угли бурые, каменные и антрациты. Классификация по генетическим и технологическим параметрам. URL: 
https://files.stroyinf-ru/Data2/1/4293772/4293772639.pdf (дата обращения: 11.01.2023). 

5 Руководство по предупреждению и тушению эндогенных пожаров на угольных шахтах Украины. КД 12.01.402-2000 от 18.12.2000. Донецк: 
Научно-исследовательский институт горноспасательного дела; 2000. 216 с. 

° Руководство по проектированию вентиляции угольных шахт. Киев: Основа; 1994. 311 c. 

7 Правила ведення гірничих робіт на пластах, схильних до газодинамічних явищ. Стандарт Мінвуглепрому України № СОУ 10.1.00174088.011- 
2005. Київ: Мінвуглепром України; 2005. 224 с. 
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Содержание минеральных примесей, судя по выходу золы из пластовых проб угля [2] и [3], может 
значительно превышать 40 96. Это, несомненно, сказывается на достоверности прогноза опасных свойств 
шахтопластов при ведении горных работ и указывает на актуальность совершенствования нормативной базы. 
Для этого необходимо разработать инженерный метод определения серы в минеральной массе углей. 

Цель представленной работы — доказать принципиальную и практически единственную возможность 
расчетного определения элементного содержания серы в минеральной массе углей для прогноза опасных свойств 
шахтопластов при ведении горных работ. Отметим, что при реализации данного подхода применяются некоторые 
стандартные показатели, зафиксированные в справочно-нормативной документации по потребительским 
качествам топлива для каждого шахтопласта. Речь идет об экспериментальных данных, касающихся содержания 
серы — общей и в горючей части топлива. Благодаря этим сведениям можно установить соотношения между 
содержанием серы в органической массе и минеральных примесях. По выходу золы можно определить 
соотношения между органической массой и минеральными примесями для обогащенных проб угля. 
Соотношение зольности пластовых и обогащенных проб дает возможность оценить содержание серы в 
удаляемых продуктах обогащения. 

Материалы и методы. Методика проведения исследований основана на тесной корреляционной связи 
(т = 0,93) содержания серы в горючей части топлива (5,) и обшей серы (SA) для углей практически всех 
шахтопластов Донецкого и Львовско-Волынского бассейнов (рис. 1). 
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Рис. 1. Зависимость содержания серы в горючей части топлива (Sr) от общей серы (5) согласно статистической обработке 
экспериментальных данных: 1 — усредняющая прямая; 
2 — биссектриса координатной сетки [2,3] 


Статистическую обработку зависимости S, от S// выполнили для 1827 шахтопластов, по которым известно 
содержание серы — общей и в горючей части топлива. 

В справочном нормативном документе [3] указаны общие сведения о потребительских качествах углей для 
734 шахтопластов, а в справочнике [2] — для 1379. 

Близость усредняющей прямой | и биссектрисы 2 координатной сетки (рис. 1) показывает, что в общем 
случае содержание серы в горючей части топлива (S-) может быть выше и ниже значения общей серы (SA). 
Разница между S/ и S. (Sj — S) в первом приближении может быть критерием соотношения элементного 
содержания серы в горючей массе и в условно негорючей части топлива. Показатель (5—5) 
для 1315 шахтопластов был отрицательным, что указывает на преимущественное содержание серы в горючей 
части массы. Для 429 шахтопластов у этого показателя были положительные значения, что свидетельствует о 
преимущественном нахождении серы в условно негорючей части топлива. Нулевая разница (5 - 5,) 
для 83 шахтопластов означала равенство содержания серы в горючей и условно негорючей частях углей для всей 
совокупности шахтопластов, рассматриваемых в статье. 

Имеется информация для 2113 шахтопластов. По некоторым из них совпадают данные в разных документах, 
и это аргумент в пользу достоверности сведений. Для четырех шахтопластов не приведены значения выхода золы 
из пластовых проб (таблица 1). По этой причине можно анализировать экспериментальные данные о выходе золы 
для 2109 шахтопластов. Из них обогащение проводилось для 1622 шахтопластов, а для 487 не проводилось. 


Техносферная безопасность 


https:// bps-journal.ru 


Безопасность техногенных и природных систем. 2024;8(1):58—72. eISSN 2541—9129 


Таблица | 


Сведения о шахтопластах, для которых не приведены данные о выходе золы из пластовых проб 


Страница 


Индекс 


Зольность, % 


Содержание серы, % 


Источник Шахта пластовой обогащенной „ са | В горючей 
источника пласта общей, S; 
пробы, А пробы, А, части, Sr 
№ 144 h = 6,3 6,3 1,9 
[3] 74 В kg = 6,1 2,8 3,0 
10-бис тз = 7,0 3,2 3,5 
[4] 106 Им. Войкова h = = = -= 


Результаты исследования. Показатели из [2] и [3] позволили рассчитать для 1827 шахтопластов содержание 


серы в условно негорючей части топлива Sy: 


S = 


H 


100- S? -(100— 4; -S, 


где А; — выход золы из обогащенных проб, 96. 
Расчетные значения S, связаны с показателем (SSi) тесной прямо пропорциональной корреляционной 


зависимостью (r = 0,91, рис. 2 a). 
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Рис. 2. Зависимость от критерия (5: — Sr) расчетных значений содержания серы ($4) по уравнению (2) в условно негорючей 


части топлива: а — для всей совокупности из 1827 шахтопластов; 6 — для 1352 шахтопластов, угли которых 
предварительно обогащались; в — для 475 шахтопластов с необогащенными углями; 1 — усредняющие прямые; 
]' — скорректированные положения усредняющих прямых; А' — условное нахождение точки для корректировки 
расположения прямой 1' и исключения получения отрицательных («диких») усредненных значений SH, х — расчетные 
значения 5н по уравнению (2) [3, 4] 


Параллельно с этим получено значительное количество отрицательных расчетных значений S, для 
802 шахтопластов. Тесная корреляция показателей и наличие «диких» (отрицательных) результатов 
свидетельствуют о возможном влиянии систематических погрешностей при определении разницы между 
содержанием серы общей (S) и в горючей части топлива (5,). Наличие отрицательных расчетных значений S, 
может указывать на то, что систематические погрешности определения критерия (5, — S.) обусловлены как 
занижением экспериментальных значений общей серы, так и завышением в горючей части топлива. 

Для раскрытия природы появления «диких» расчетных значений содержания серы в условно негорючей части 
топлива по уравнению (2) рассмотрели методики определения показателей, входящих в эту зависимость. 

Известны несколько подходов к определению общей серы в углях — от классического гравиметрического 
метода Эшка до современных видов инструментального анализа [1]. Метод Эшка, предложенный в 1874 году, 
все еще широко используется в мировой практике. Он регламентируется ГОСТом. Навеску угля сжигают путем 
спекания со смесью оксида магния и безводного углекислого натрия (соотношение 2:1 по массе) в окислительной 
атмосфере при температуре (800 + 250) °С. Органическая масса топлива сгорает, и сера превращается в сульфаты 


натрия и магния. Их переводят в раствор. Сульфат-ион количественно осаждают в солянокислой среде 
хлористым барием в виде сернокислого бария Ваз Ол. По массе Ва5 О. рассчитывают общее содержание серы. 
Недостатки этого подхода: большие затраты времени (не менее 6-7 часов) и возможность потери части 
органической серы при спекании со смесью [1]. В большинстве случаев метод применяется при возникновении 
разногласий как арбитражный. Общую серу в стандартных образцах топлива, предназначенных для градуировки 
анализаторов, определяют только методом Эшка. 

Для практических целей рекомендуется регламентированный ГОСТом? ускоренный стандартный метод 
определения общей серы. Навеска топлива сжигается в потоке кислорода или воздуха в трубчатой печи при 
1350 °С. Образуются зольный остаток и газообразные продукты сжигания, в состав которых входят оксиды серы 
(в основном 502) и хлора. Анализ выполняется за 40-60 минут. Этот метод использовали при определении 


$ Правила ведення гірничих робіт на пластах, схильних до газодинамічних явищ. Стандарт Мінвуглепрому України СОУ 10.1.00174088.011- 
2005. Київ: Мінвуглепром України; 2005. 224 с. 

? ГОСТ 2059-95 (ИСО 351-96). Топливо твердое минеральное. Метод определения общей серы сжиганием при высокой температуре. Москва: 
Стандартинформ; 2008. 14 с. URL: https://files.stroyinf.ru/Data2/1/4294832/4294832840.pdf (дата обращения: 11.11.2023). 
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содержания общей серы для установления потребительских качеств углей, экспериментальные значения 
приведены в справочно-нормативных документах [2,3]. Ускоренный метод определения общей серы 
рекомендуется для контроля за содержанием серы в малозольных топливах. Его нежелательно использовать при 
анализе высокозольных топлив со значительным содержанием соединений кальция. Из-за разнообразия состава 
минеральной массы и соединений серы нельзя быть полностью уверенным в том, что часть серы не останется в 
золе в виде CaSO4. Как следствие, результаты получаются заниженными [1]. 

Для большинства шахтопластов после обогащения угля выход золы составляет менее 10 %, но этот показатель 
бывает и значительно выше (рис. 3). 
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Рис. 3. Зависимость выхода золы из обогащенных проб угля (Ao°) от выхода золы из пластовых проб (An°) согласно 
экспериментальным данным [2, 3] 


По этой причине даже после обогащения значительной части проб нет оснований считать все угли малозольными. 
Для 487 шахтопластов, согласно данным [2, 3], обогащение не проводилось и выход золы во многих случаях превышал 
10 96, что также не позволяет относить угли этой совокупности шахтопластов к малозольным. 

Соединения кальция есть в минеральных примесях углей всех шахтопластов. Об этом свидетельствует 
типизация состава зол (таблица 2), согласно которой содержание CaO в общем случае может достигать 40 %. 


Таблица 2 
Типизация состава зол углей по данным [2] и [3] 
2 Пределы колебания состава, % 
Типы зол углей : 

АБО; SiO» Fe303 CaO 

Кремнистые 8—30 40—70 до 20 до 20 

Глиноземные 30-45 40-55 до 20 до 20 

Железистые 10—20 30—55 >20 до 20 

Известковые 5—20 15—40 5—20 20—40 
Одна из причин получения «диких» расчетных значений 5, по уравнению (2) — занижение S^ при 


определении ускоренным методом. Такой вывод обосновывается: 

— анализом выхода золы из обогащенных и пластовых проб; 

— возможным содержанием соединений кальция во всех типах зол. 

Второй фактор получения отрицательных расчетных значений S, — завышение доли горючей части топлива 
с использованием в качестве множителя (100—A,°) в уравнении (2). Содержание общей серы 5/ относится как к 
органической части, так и к минеральной массе (ММ) топлива. По этой причине более логичным было бы 
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использовать в уравнении (2) вместо выхода золы А, значений ММ, которые для рассматриваемых шахтопластов 
экспериментально не определялись. Согласно известным эмпирическим зависимостям, значения ММ всегда 
превышают выход золы на сухое состояние. Учет такого превышения требует дополнительных (недостающих) 
сведений о гидратной влаге, выделяющемся диоксиде углерода при полном разложении карбонатов, пиритной и 
сульфатной сере. Для приблизительной оценки содержания минеральной массы иногда используют соотношение 
MM = 1,1 4 [1]. 

Анализ показал причины систематического занижения произведения 100-S/ и завышения разности 
(100—4,*) S. при определении расчетных значений S, по уравнению (2). 

Наличие отрицательных значений S, для совокупности из 802 шахтопластов (рис. 2 а) дает возможность 


оценить среднее занижение критерия в = КА ‚ при котором получены «дикие» результаты: 


п 


2.45. -5,) 
Fy sme. 


0,42, % (3) 


Здесь n = 802 — количество шахтопластов, для которых получены отрицательные значения Su; z(s? =S; ) = 
А 


сумма отрицательных значений 5; — S. для совокупности из и = 802 шахтопластов. 


Среднему отрицательному значению ig — S, = – 0,42% соответствует среднее отрицательное значение Sa 
n 
25. 1682,18 
S 5t = ——=—2,10,%, (4) 
п 802 
n 
rye > 5; — Сумма отрицательных расчетных значений S, для совокупности из п = 802 шахтопластов. 


ГА 


Чтобы исключить получение отрицательных усредненных значений Su, необходимо положение прямой 1 


(рис. 2 а) скорректировать с учетом координат 57 — S. = —0,42 % и Sa = —2,10 %. Скорректированная прямая |' 


должна проходить не ниже точкиА' с координатами A = =-0,42% и S,=+2,10%. Прямо 


пропорциональная зависимость $. OT (5—5) (r=0,91) указывает на TO, что истинное положение 
скорректированной прямой 1' будет близко к ее параллельности с прямой 1. Использование в инженерных 
расчетах значений S,, определенных по уравнению 1' (рис. 2 а), гарантирует, как минимум, отсутствие их 
отрицательных усредненных значений. 

Эмпирическое уравнение 1' получено на основании экспериментальных данных для 1827 шахтопластов. Из 
них для проб углей 1352 шахтопластов проводилось предварительное обогащение. Выход золы из пластовых 
проб углей остальных 475 шахтопластов был несколько ниже, поэтому для них обогащение не 
предусматривалось. 

Предварительное обогащение существенно повлияло на выход золы по сравнению с необогащенными 
пробами (рис. 3). Следовало выяснить, как точность расчетного определения серы по уравнению (2) в условно 
негорючей части топлива зависит от обогащения. Для этого рассмотрели отдельные характерные совокупности 
шахтопластов. К одной совокупности из 1352 вариант отнесли шахтопласты, угли которых подвергались 
предварительному обогащению (рис. 2 6). До обогащения выход золы из пластовых проб составлял 11,0-48,2 %, 
а после — 3,7—16,7 % (рис. 3). 

Во второй совокупности из 475 вариант рассмотрели шахтопласты, угли которых предварительно не 
обогащались (рис. 2 в), так как условно их относили к малозольным. Выход же золы из необогащенных углей 
этой совокупности шахтопластов колебался от 2,4 % до 28,1 % (рис. 3). 

Выход золы для обеих рассматриваемых совокупностей шахтопластов во многих случаях значительно 
превышал 10%. Диапазоны выхода составляли 3,7—16,7 % и 2,4—28,1 % соответственно. Это указывает на 
отсутствие в ряде случаев для обеих совокупностей обоснованного применения альтернативного ускоренного 
стандартного метода определения общей серы 55, [9]. В результате по большинству рассматриваемых 
шахтопластов расчетные значения 5/ оказались заниженными. Вследствие этого получены и «дикие», 
рассчитанные по уравнению (2), значения Su соответственно для 558 и 244 шахтопластов обеих (2-й и 3-й) 
совокупностей (таблица 3). 
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Средние отрицательные значения 5? —5 для шахтопластов 2-й и 3-й совокупностей были рассчитаны по 


уравнению (3). Они равнялись соответственно —0,42 % и —0,41 %. Это указывает на практически одинаковые 
средние погрешности занижения расчетного содержания общей серы альтернативным ускоренным методом для 
обеих совокупностей при расчете $ по уравнению (2). По указанным обстоятельствам ускоренный метод 
определения обшей серы и возможные занижения ее истинных значений не относятся к различиям между 2-й и 
3-й совокупностями шахтопластов. 

Логика уравнения (2) показывает, что разные расчетные значения Sy 2-й и 3-й совокупностей можно получить 
только с учетом разного выхода золы из обогащенных и пластовых проб углей. Для удобства анализа 
потенциального влияния выхода золы из обогащенных и необогащенных проб углей уравнение (2) представили 
в виде слагаемых. 

Для обогащенных проб (совокупность 2): 


1 1 
se =P Sf TUS ess (5) 


Для необогащенных пластовых проб (совокупность 3): 


u 1 1 п п 
S, N a онота (6) 
AS AS I I 


п II 


Сравнивая слагаемые уравнений (5) и (6), определили направленность их влияния на расчетные значения Sy. 
d d 
Первое слагаемое (100/ Ае ) .5, уравнения (5) заведомо больше аналогичного слагаемого (100/ A ):5, 


уравнения (6). Это определили по бесспорному превышению выхода золы из пластовых проб по сравнению с 
выходом из обогащенных (Án? > Ao’). Неравенство позволяет утверждать, что значения общей серы, определенные 
ускоренным альтернативным методом, будут в большей степени занижаться для пластовых проб. По этой 
причине значения 5; уравнения (5), как минимум, не меньше, чем в (6). 


Вторые слагаемые со знаком минус уравнений (5) и (6) ((100/ 4€). 5? и (100/ 4;).57) характеризуют 


увеличение доли содержания серы в горючей части топлива после обогащения углей по сравнению с 
содержанием в пластовых пробах до обогащения. Под воздействием процессов обогащения выход золы 
сокращается в несколько раз (4, < Аи‘, рис. 3). В результате увеличивается доля горючей части топлива 


обогащенных углей (100/42 ) по сравнению с ee долей в пластовых пробах (100/4° ). Это, a также 


устойчивость соединений серы обусловливают рост относительного (процентного) содержания серы S; В 
горючей части топлива после обогащения. Эти значения всегда превышают процентное содержание серы в 
пластовых необогащенных пробах, т.е. соблюдается неравенство S. Sp Отрицательные слагаемые 


уравнений (5) и (6) ((109/ 45). 5? и (100/4;)-s") при любом их соотношении составляют лишь некоторую 


часть от содержания серы в горючей массе. Вследствие этого суммарные значения второго и третьего 
слагаемых уравнений (5) и (6) будут положительными. 


Рассмотрим неравенство слагаемых B (5) и (6): 5° > SF u (100/ 4; )- s? и (100/4 )-5/. Это говорит о том, 


что обогащение ведет к процентному увеличению содержания серы B негорючей части топлива (Sy) и ee 
тесной, почти функциональной корреляционной зависимости (ғ = 0,94) от критерия 55, (рис. 2 б). Сильная 
корреляция Su и SÉ — S. для обогащенных углей шахтопластов 2-й совокупности обусловлена значительным 
принудительным сокращением содержания условно негорючих (минеральных) примесей. Об этом говорит 
сокращение диапазона выхода золы из пластовых необогащенных проб с 11,0-48,2 % до 3,7—16,7 % после 
обогащения. Для рассматриваемой совокупности шахтопластов в большинстве случаев (около 90 %) выход 
золы составлял менее 10 %, что указывает на незначительную долю содержания минеральных примесей в 
углях после обогащения. Обогащением в значительной степени объясняется приближение массовой доли 
горючей части топлива к органической массе. Корреляционная зависимость 1 (рис. 2 б) относится к 
незначительной доле (менее 10 %) сохранившихся после обогащения минеральных (условно негорючих) 
примесей. Относительно высокий коэффициент регрессии уравнения 1 (13,15) указывает на диспропорцию 
показателей: 

— сокращения выхода золы после обогащения углей; 

— содержания серы в сохранившейся части условно негорючих (минеральных) примесей [10]. 
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Как видно из уравнения (5), высокая корреляционная зависимость 5, от показателя Si — 5, для шахтопластов 
2-й совокупности обеспечивалась относительным постоянством выхода золы Ао° для большинства шахтопластов 
в пределах 3,7—10,0 %. Это фактически исключает существенное влияние непредсказуемых колебаний доли 
условно негорючих (минеральных) примесей в более широких диапазонах. Значения Sy? для 2-й совокупности 
шахтопластов, рассчитанные по уравнению (5), определяются экспериментальными данными 5/ и S. Это 
подтвердили результаты статистической обработки (рис. 2 6). Влияние колебаний выхода золы Ао’ после 
обогащения углей свелось до минимума, и расчетные значения Su? в большей степени зависели от двух 
показателей — S/ и S,°. 

Такой устойчивой связи между 5, и (S-S) не наблюдается при отсутствии процессов 
обогащения (рис. 2 в). Это выразилось в более низком коэффициенте корреляции (r = 0,56) и снижении 
коэффициентов регрессии эмпирических уравнений | с 13,15 (рис. 2 6) до 3,50 (рис. 2 в). 

При расчете значений 5." по уравнению (6) для 3-й совокупности шахтопластов диапазон 
непредсказуемых колебаний значений выхода золы 4: расширился до 2,4 %-28,1 %. Это значительно 
больше, чем при стабильном выходе золы из обогащенных проб углей 2-й совокупности. Отметим усиление 
влияния дополнительного фактора, связанного с существенными колебаниями выхода золы из пластовых 
проб в более широких интервалах. Это привело к ослаблению корреляционной зависимости | до значения 
r = 0,56 (рис. 2 в). 

Подавляющее большинство углей 3-й совокупности шахтопластов по критерию 54—-5,<0 должно 
относиться к малозольным. Из 475 шахтопластов всей рассматриваемой совокупности для 454 (т. е. для 95,6 %) 
значения критерия были отрицательными. По отрицательным значениям критерия 5; — 5, основное процентное 
содержание серы необходимо было условно отнести к горючей массе углей этих шахтопластов. Определение 
Sa согласно уравнению (6) для 244 шахтопластов дало «дикие» значения. Для остальных 231 шахтопласта 
расчетные значения S, были положительными (210 вариант) или равнялись нулю (21 варианта). Исследование 
выявило примерно равное количество околонулевых расчетных значений Sy. Это подтверждает наличие 
систематических погрешностей экспериментально определенных параметров, входящих в (6). 
Доказательством служит также примерное равенство по абсолютной величине средних отрицательных 
(0,71 %) и средних положительных (+0,80 %) расчетных значений Sy. Это дает возможность обоснованно 
вносить количественные корректировки в расчеты, чтобы исключить получение «диких» результатов. 

Сравнение расчетных значений содержания серы в негорючей части топлива для 2-й и 3-й совокупностей 
демонстрирует значительное отличие исходных параметров показателей, которые использовались при 
определении Sy из," по уравнениям (5) и (6) соответственно. Основная разница обусловлена 
предварительным обогащением углей 2-й совокупности. Обогащение непропорционально изменило 
соотношение выхода золы и содержания серы в горючей части топлива. Исходные соотношения 
формировались в природных условиях и связаны с индивидуальной геологической историей образования 
каждого шахтопласта, а именно: 

— последовательным накоплением осадков; 

— погружением и преобразованием исходного вещества на предшествующих метаморфизму стадиях 
торфообразования и диагенеза. 

В углях 3-й совокупности шахтопластов сохранились исходные сложившиеся природные соотношения между 
выходом золы из пластовых проб, содержанием общей серы и некоторой ее части в горючей массе топлива. 

Обогащение углей 2-й совокупности шахтопластов и сокращение выхода золы (в большинстве случаев 
ниже 10%) предопределили искусственную близость всей совокупности по содержанию минеральных 
(условно негорючих) примесей. В результате фактического устранения влияния выхода золы из 
обогащенных проб для большинства шахтопластов по уравнению (5) получены расчетные значения S", 
которые в значительной степени непосредственно связаны с критерием (Sf — Sr). Это и предопределило их 
высокую корреляционную взаимозависимость (рис. 2 6). 

Главное отличие углей 2-й и 3-й совокупностей шахтопластов — значения выхода золы из обогащенных 
и пластовых проб углей. Для детального представления о влиянии выхода золы из обогащенных и 
необогащенных проб углей на расчетные значения 5, дополнительно рассмотрели четыре совокупности 
шахтопластов с углями разной степени глубины обогащения (рис. 4). 
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Рис. 4. Зависимость расчетных значений содержания серы (Su) в условно негорючей части топлива от критерия (5 — Sr) [2, 3] 


По выходу золы из обогащенных проб шахтопласты разделили на две группы: менее 10 % и более 10% 
(рис. 3). Рассмотрели соответствующие группам совокупности — 4-ю (рис. 4 а) и 5-ю (рис. 4 б). 

Большинство шахтопластов (1181 % или 87,4 %) по выходу золы менее 10 % из обогащенных проб отнесли к 
4-й совокупности. В 5-й совокупности шахтопластов с выходом золы более 10% рассмотрели расчетные 
значения Sy для 171 шахтопласта, что составило 12,6 % от общего количества шахтопластов (1352), угли которых 
обогащались. 

Для обеих совокупностей шахтопластов после обогащения получены тесные прямо пропорциональные 
корреляционные зависимости расчетных значений S, от критерия 5/ — S. Коэффициенты корреляции — 0,95 
и 0,92 соответственно. Коэффициенты регрессии этих зависимостей (рис. 4а, б) равнялись 13,84 и 7,50 
соответственно. Почти двукратная разница между эмпирическими коэффициентами регрессии свидетельствует 
о существенных различиях между углями 4-Й и 5-й совокупностей. Ее можно объяснить: 

— количественным и качественным составом минеральных примесей в пластовых пробах; 

— разной эффективностью обогащения. 

Влияние указанных факторов требует отдельного детального исследования [11]. 

Аналогичные совокупности шахтопластов сформировали по выходу золы из пластовых (необогащенных) 
проб менее и более 10 %. Для 6-й совокупности из 149 шахтопластов установили зависимость Sy от критерия 
5—5; при выходе золы из необогащенных (пластовых) проб менее 10 % (рис. 4 в). В 7-й совокупности из 326 
шахтопластов рассмотрели аналогичные зависимости расчетных значений $, от критерия (S/— Sr) при выходе 
золы более 10 % (рис. 4 2). 

Полученные прямо пропорциональные эмпирические зависимости для 6-й и 7-й совокупностей шахтопластов 
с разным выходом золы из необогащенных (пластовых) проб значительно различаются. Отличия обусловлены 
как теснотой корреляционных связей, так и коэффициентами регрессии эмпирических уравнений. 
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Выход золы из необогащенных (пластовых) проб менее 10 % указывает на низкую зольность углей 6-й 
совокупности шахтопластов. В этом случае корреляционная связь довольно высока (r = 0,79). Для более зольных 
углей 7-й совокупности шахтопластов корреляция значительно слабее (7 = 0,56). 

Более чем в три раза различаются коэффициенты регрессии эмпирических уравнений 6-Й и 7-й совокупности 
шахтопластов. Они соответственно равны 8,27 (рис. 4 в) и 2,53 (рис. 4 2). 

Выявленные существенные отличия между эмпирическими зависимостями 6-й и 7-й совокупностей шахтопластов 
по тесноте корреляционных связей и коэффициентов регрессии указывают на разные особенности условий 
накопления серы в минеральной (условно негорючей) части топлива для углей с разным уровнем выхода золы. 

Низкий уровень выхода золы (менее 10 %) свидетельствует о незначительном содержании минеральных 
примесей. Между расчетным содержанием серы в условно негорючей части и критерием (S/ — Sr) (рис. 4 в) 
фиксируется довольно высокий коэффициент корреляции (г= 0,79). Это значит, что сера формировалась 
параллельно в минеральных примесях и в органической массе на разных стадиях образования угля. Отсутствие 
такой тесной корреляции для высокозольных углей указывает на разные условия формирования содержания серы 
в минеральных примесях и в органической массе 7-й совокупности шахтопластов (рис. 4 г). 

Обсуждение и заключение. Согласно эмпирическим уравнениям (рис. 4) попарно рассмотрели совокупности 
шахтопластов, которые при прочих равных условиях отличались по двум факторам (таблица 3): 

— обогащались или не обогащались угли; 

— выход золы из обогащенных или необогащенных углей был менее или более 10 %. 

Обогащение существенно меняет соотношения между: 

— органической массой и минеральными примесями; 

— горючими и условно негорючими частями топлива. 

По этой причине 1-я, 2-я, 4-я и 5-я совокупности шахтопластов и соответствующие им эмпирические 
уравнения не отражают истинное природное содержание серы в негорючих частях ископаемых углей. Их нельзя 
использовать для прогноза опасных свойств шахтопластов при ведении горных работ. 

Содержание серы в условно негорючей массе необогащенных проб для 3-й и 7-й совокупностей занижается 
из-за наличия высокого выхода золы и применения альтернативного стандартного метода определения 
содержания общей серы. Область применения соответствующих этим совокупностям эмпирических уравнений 
ограничена шахтопластами с низкозольными углями и незначительным содержанием соединений кальция в 
минеральных примесях. 

Наиболее достоверные результаты определения расчетного содержания серы в условно негорючей массе 
получены для 6-й совокупности шахтопластов. Практическое применение соответствующего эмпирического 
уравнения рекомендуется при прогнозе опасных свойств шахтопластов с выходом золы менее 10 %. 

Можно улучшить предложенный метод определения серы в минеральных примесях и расширить область его 
применения. Для этого необходимо установить соотношения: 

— содержания серы в органической и горючей массе углей; 

— минеральных примесей и выхода золы (с определением ее химического состава). 

Результаты таких исследований можно использовать для совершенствования нормативов по безопасному 
ведению горных работ с учетом содержания серы как в органической массе, так и в минеральных примесях 
ископаемых углей. 
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Аннотация 


Введение. Вопросы обеспечения оптимальных режимов расхода топлива являются предметом постоянного 
исследования в области эксплуатации автотранспортных средств. Работы в области снижения расхода топлива 
в большинстве случаев посвящены исследованию конструкции силовых установок, трансмиссии или 
аэродинамических особенностей кузова автомобилей. Вместе с тем вопросы определения оптимальных законов 
управления неустановившимся движением с точки зрения синтеза алгоритмов управления движением 
автомобиля не рассматривались. Для оптимизации передаточного числа автомобилей при неустановившимся 
движении авторами работы проведены аналитические расчеты, позволяющие упростить синтез алгоритмов 
управления движением, что позволяет снизить расход топлива. Целью работы являлось определение 
оптимального передаточного числа, позволяющего решить задачу синтеза управления транспортным средством 
для снижения расхода топлива. 

Материалы и методы. Рассматриваемый в работе синтез алгоритмов управления движением автомобиля 
основан на применении к инвариантным признакам действительного движения игольчатой вариации 
Л.С. Понтрягина. Использован аналитический метод оценки энергетической эффективности эксплуатационных 
характеристик автомобиля, в основе которого лежит определение оптимального передаточного числа 
трансмиссии автотранспортных средств с учетом наименьшего расхода топлива. Представленный метод 
учитывает величину крутящего момента, передаваемого от коленчатого вала двигателя на элементы трансмиссии, 
который, в зависимости от мощности двигателя, реализуется в виде силы тяги на колесах автомобиля. 
Результаты исследования. Построен закон оптимального изменения передаточного числа трансмиссии при 
разгоне автомобиля за минимальное время. Решена задача определения оптимального передаточного числа 
трансмиссии автомобиля в случае движения автомобиля с постоянной скоростью и постоянной подачей 
топлива, а также при разгоне автомобиля до заданной скорости при постоянной подаче топлива, когда 
соблюдается условие = = const. Результатом рассмотренного случая применения оптимального закона 
изменения передаточного числа трансмиссии является минимизация расхода топлива при ограничениях на 
ускорение (силу тяги) и скорость движения автомобиля. 

Обсуждение и заключение. Применение авторами аналитического метода для определения передаточного числа 
трансмиссии автомобиля, а также использование данного метода в практических расчетах для автомобиля с 
заданными характеристиками показывает возможность решения задачи синтеза управления транспортным 
средством с применением математического аппарата. Это подтверждается построенной графической 
зависимостью по результатам проведенных расчётов. Рассмотренные случаи движения позволили определить 
аналитические зависимости оптимального передаточного числа и скорости автомобиля. Исходные данные, 
полученные по аналитическим зависимостям, применимы для автомобилей со смешанным режимом управления. 


Ключевые слова: передаточное число трансмиссии, расход топлива, энергетическая эффективность 


Благодарности. Авторы выражают благодарность редакции и рецензентам за внимательное отношение к 
статье и указанные замечания, устранение которых позволило повысить ее качество. 


© Косенко E.E., Марченко Ю.В., Марченко Э.В., Попов C.H., Изюмов A.H., 2024 


Машиностроение 


https:// bps-journal.ru 


Косенко Е.Е. и др. Оптимизация передаточного числа трансмиссии по расходу топлива 


Для цитирования. Косенко E.E., Марченко Ю.В., Марченко Э.В., Попов C.M., Изюмов А.И. Оптимизация 
передаточного числа трансмиссии по расходу топлива. Безопасность техногенных и природных систем. 
2024;8(1):73—81. https://doi.org/10.23947/2541-9129-2024-8- 1-73-81 


Original article 


Optimization of the Transmission Ratio by Fuel Consumption 


Evgeny E. Kosenko®, Julianna V. Marchenko®, Edward V. Marchenko®, Sergey I. Popov®, 
Andrey I. Izyumov® 


Don State Technical University, Rostov-on-Don, Russian Federation 


al23lok@mail.ru 


Abstract 

Introduction. The issues of ensuring optimal fuel consumption modes are the subject of constant research in the field of 
vehicle operation. Works in the field of reducing fuel consumption in most cases are devoted to the study of the design 
of power units, transmission or aerodynamic features of the body of cars. At the same time, the issues of determining 
the optimal laws of control of unsteady movement from the point of view of the synthesis of algorithms for controlling 
the movement of a car have not been considered. To optimize the transmission ratio of cars with unsteady movement, 
the authors of the work carried out analytical calculations to simplify the synthesis of motion control algorithms, which 
in turn allowed reducing fuel consumption. The aim of the work was to determine the optimal transmission ratio, which 


allowed solving the problem of synthesizing vehicle control to reduce fuel consumption. 

Materials and Methods. The synthesis of algorithms for controlling the movement of a car considered in the work was 
based on the application of the needle variation of L.S. Pontryagin to invariant features of the real movement. An 
analytical method was used to estimate energy efficiency of vehicle performance, which was based on determination of 
optimal transmission ratio of motor vehicles taking into account minimum fuel consumption. The presented method 
took into account the amount of torque transmitted from the engine crankshaft to the transmission elements, which, 
depending on the engine power, was realized in the form of traction force on the wheels of the car. 

Results. The law of optimal change in the transmission gear ratio during acceleration of the car in a minimum time was 
built. The problem of determining the optimal transmission ratio of the vehicle in the case of driving the vehicle at a 
constant speed and constant fuel supply and in the case of accelerating the vehicle to a given speed at a constant fuel 
supply, when the condition = = const was met, was solved. The result of the considered case of applying the optimal law 
of change in the transmission gear ratio was the minimization of fuel consumption under restrictions on acceleration 
(traction force) and speed of the car. 

Discussion and Conclusion. The use by the authors of the analytical method for determining the transmission ratio of a 
car, as well as the use of this method in practical calculations for a car with given characteristics, showed the possibility 
of solving the problem of synthesizing vehicle control using a mathematical apparatus. This was confirmed by the built 
graphical dependence based on the results of the calculations. The considered cases of movement made it possible to 
determine the analytical dependencies of the optimal transmission ratio and the speed of the car. The initial data 
obtained by analytical relationships are applicable for cars with a mixed control mode. 
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Введение. При движении автомобильного транспорта, с учетом его перемещения в условиях городских 
застроек, режим управления транспортным средством является смешанным с чередованием разгона и 
торможения, а также с необходимостью различного маневрирования. Подобное управление приводит к 
повышенному расходу топлива. Уменьшение расхода топлива в Таких условиях, как правило, связано с 
необходимостью увеличения времени совершения маневров [1]. Это вызывает необходимость поиска 
оптимальных законов управления неустановившимся движением автомобиля. В работе [2] энергетические 
характеристики перемещения транспортного средства определяются в первую очередь величиной 
передаваемого от коленчатого вала двигателя крутящего момента и реализуемой на движителях силы тяги, 
которые зависят от мощности двигателя и значения передаточных чисел. Необходимость определения 
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оптимальных значений передаточных чисел трансмиссии транспортного средства обусловлена уменьшением 
значений топливных характеристик [3], а также уменьшением величины энергии, которая тратится при 
маневрировании автомобиля в случае неустановившегося движения. Передаточные числа определяются в 
соответствии с заданной характеристикой двигателя и в зависимости от скорости перемещения транспортного 
средства, имеющей переменные значения [4]. 

При рассмотрении применяемых подходов для решения задач, связанных с оптимизацией, следует сказать 
об их многообразии. Большая часть таких подходов использует методы оптимального управления и 
вариационного исчисления. Авторы в [5] использовали принцип максимума Л.С. Понтрягина с определением 
аналитических зависимостей для гибридных автомобилей. Подобный подход рассматривается и в работе [6], на 
основе которого предлагается построение закона оптимального изменения передаточного числа трансмиссии 
при разгоне автомобиля за минимальное время, а также определение топливной экономичности в ездовом 
цикле [7]. Анализ проведенных исследований показал значительный теоретический задел в области 
определения оптимального передаточного числа при различных режимах эксплуатации автомобиля. Вместе с 
тем, следует сказать, что оптимальный закон изменения передаточного числа трансмиссии должен 
определяться с учетом минимизации расхода топлива при ограничениях на ускорение (силу тяги) и скорость 
движения автомобиля. Авторами представленной работы предлагается для определения оптимальных законов 
управления неустановившимся движением (с точки зрения синтеза алгоритмов управления) применить метод 
оптимизации, в основе которого лежит применение к инвариантным признакам действительного движения 
игольчатой вариации Л.С. Понтрягина. Предлагаемый подход авторами ранее проведенных исследований не 
рассматривался. 

Целью работы является определение оптимального передаточного числа, позволяющего решить задачу 
синтеза управления транспортным средством для оптимизации расхода топлива. 

Материалы и методы. При расчете мгновенного расхода топлива автомобиля используем следующую 
зависимость: 

Ст - quK, Ky N,, (1) 
где фу — показатель удельного расхода топлива, соответствующий максимальному значению № = Nmax, 
К, — значение коэффициента, определяемого с учетом скоростного режима работы двигателя; Ky — значение 
коэффициента, определяемого с учетом степени загрузки двигателя; № — мощность двигателя в текущий 
момент времени. 

Коэффициенты К, и Км определяли с учетом эмпирических зависимостей, К» — с учетом скоростного 
режима работы двигателя. Определим его по формуле: 
2 


o’ : , (2) 
On On 


где а» = 1,27; b, = 0,94; c, = 0,67; о, — текущее значение угловой скорости коленчатого вала; оу — угловая 
скорость вращения коленчатого вала, соответствующая максимальной мощности двигателя Nmax- 
Величина о, связана с передаточным числом трансмиссии i и скоростью движения автомобиля X 


зависимостью: 
0, = (i: X)/ r,. (3) 
Тогда формула (2) примет вид: 
b | 6 РЕ 
Ko = lo — ({-х)+——®—(1-%).. (4) 
Onl ( Oyr, ) 


Коэффициент Ky зависит от степени загрузки двигателя 5. Значение = определяли отношением мощности 
двигателя в текущий момент времени независимо от режима движения транспортного средства. Принимали В 
расчет: текущее значение угловой скорости коленчатого вала We, мощность двигателя Nec в момент полной 
подачи топлива и такое же значение угловой скорости коленчатого вала We: 


s=—, (5) 


Значение текущей мощности двигателя № определяли дифференциальным уравнением движения 
автомобиля по оси x, имеющего вид [5]: 


06. N : 
Sg- дейт y-k. (6) 
g mgx 

Здесь б —— значение коэффициента, определяемого с учетом затрат мощности с учетом приращения 


кинетической энергии вращающихся масс двигателя, трансмиссии и колес при разгоне автомобиля; 
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"nr — коэффициент полезного действия трансмиссии; т — полная масса автомобиля; 
V = Ѓсоѕа + sina — коэффициент сопротивления движению; К — коэффициент обтекаемости, определяющий 
силу сопротивления воздуха. 

Коэффициент 8 в представленной работе определили по следующей зависимости à = 1 + үі?, в которой 
значение у зависит от конструктивных параметров автомобиля и изменяется в пределах у = 0,001—0,003. Тогда 
из уравнения (6) можно получить выражение для текущего значения мощности двигателя: 

_ тах | 1+& : 

Nr & 

Значение мощности двигателя Nec в момент максимальной подачи топливной смеси определили с помощью 
известной зависимости (формула Лейдермана) для внешней характеристики двигателя, имеющей вид: 

2 
B. dispen 2| |, (8) 
N On On 


N 


e 


ž+y(t)+ k? |. (7) 


[0) 
=N, max 


N 


ec 


Эмпирические коэффициенты A, В, С, зависящие OT параметров двигателя, могут определяться, например, 
по зависимостям, приведенным в работе [5]. После преобразований с учетом (3) выражение (8) принимает вид 


Мы = ЕЁ(ах+Ы?х? ci? X ), (9) 


где коэффициенты а, Ь и с определили по формулам: 


N B C 
Е". а= Д: b= -c-— 
Onl, бла Onl, = (уту 
Эмпирическая зависимость для коэффициента Ky имеет вид: 
Ков ее ody й (10) 
где константы ау = 3,27; by = — 8,22; cy = — 9,13; dy = 3,18 для бензиновых двигателей; ay = 1,20; by = 0,14; 
см = 1,80; dy = — 1,46 для дизельных двигателей. 


Степень загрузки двигателя = в зависимости (10) с учетом (5), (7) и (9) определили по формуле: 
2 

mgaoyr, | 1+ .. : 

ЕМА CITY ea y(t) +i? 

N maT 8 

г "au PR . (11) 
ai + bi"x-cix 

В качестве целевой функции использовали величину расхода топлива за время fk, в течение которого 


осуществляется движение автомобиля на траектории хү) < х, C переменной скоростью x= f(t). Задача 


оптимизации решали в соответствии со следующим алгоритмом: устанавливали закон, в соответствии с 
которым определяется изменение передаточного числа трансмиссии i( x ) ‚ И соответствующую ему траекторию 


x(t) C таким учетом, чтобы целевой функционал принимал минимальное значение: 


tk 
Ј = [оа > min, (12) 
0 
передаточное число, текущие скорость и мощность удовлетворяли ограничениям: 
Ї> igs XS Хау, М, < N naxs (13) 


где іо — передаточное число главной передачи. 
Результаты исследования. Результаты проведенных расчетов были использованы для случая 
непрерывного изменения передаточного числа трансмиссии. Тогда, пренебрегая инерцией вращающихся масс 


(при y = 0), оптимальное значение її, если оно существует, находится из условия: 
(Gri) (ок, оку К i 
D Os солу ш (14) 
ді ði а Ky ді 


Полученное выражение (15) в качестве результатов исследования применили для двух предельных случаев. 


Первый случай — движение автомобиля с постоянной скоростью X = const и постоянной подачей топлива, 
когда € и Ky = const. Поэтому второе слагаемое в уравнении (14) равно нулю, что позволяет получить простую 
формулу для определения оптимального передаточного числа трансмиссии: 

P by 7, K 


i = 15 
i 2с„Х (13) 
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При минимизации времени разгона до заданной скорости оптимальное передаточное число трансмиссии 
определили по формуле [5]: 


Bo yr, 3AC 
iP, = Мк |+ (16) 
LEE 7 В? 
Второй случай — разгон автомобиля ДО заданной скорости при постоянной подаче топлива, когда 


соблюдается условие є = const. В этом случае значение коэффициента Ky = const и оптимальное передаточное 
число трансмиссии определили по формуле из решения уравнения: 


OK, | 
+ К, Е 0. 
д Oi 
После преобразований с учетом выражения (4) для К„ получим: 
b b | 
ж ig? | as (iX)  — M (xy = (17) 
One ( yr, ) Onl, (Gyr, ) д 
ox 


Производную a в уравнении (17) определили из условия € = const c использованием выражения (11). 
i 


При этом будем пренебрегать, Kak и ранее, инерцией вращающихся масс (y = 0) и силой сопротивления воздуха 
(кг = 0) . Тогда связь между скоростью движения автомобиля и передаточным числом трансмиссии установим 
(при € = const) из решения следующего уравнения: 
За „з .. z(t 
c? x -bi?x-ai (t) 
ed 


== 0, 
где введены обозначения: 


" ; 
= тт „уу. чу (и |. 
т89 yr; 8 


Полученное уравнение относительно скорости имеет решение, если его дискриминант больше или равен нулю. 


3 
Aci 


ed 
Тогда для передаточного числа трансмиссии должно выполняться условие: 


D = 21 + Дасі“ — 


z(t) = 0. 


oe 18 
"S (B /4C +A) = 


Соответствующее этому условию значение скорости определили по формуле: 


При разгоне автомобиля предельное (минимальное) значение передаточного числа трансмиссии Тит 


ox 
определили из условия (18), что позволяет получить достаточно простое выражение для производной ài H 
1 


распространить его на всю траекторию движения: 
дї b B 
a we бй ТЕ (20) 


тіп тіп 


Подставляя полученные выражения в уравнение (17), после преобразований получим уравнение для 
оптимального значения передаточного числа трансмиссии: 


bo 2 бы 1-Я? |+ | a, % (jg). Ч z(i-xy [A =0, 
Oy, ( Oy, ) Oy. (yr, ) ôi 


Bb, cB? Boyn b, B (2 | В 
© B+ a 2 Fc Oy, =0, 
с" 2С? В 2Ci. p © 2С В [0] а? 2СР N'k 


min 


b, 2сВ | Фут, | 2Са Б, teal B ) Oy 0, 
is is В В 2С, 


С 
RU D ——- 
L C p e Ww 
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Ч Ща 
- 2 
2с„В. Ri. 
—b + оү — 
С Р 


тт 


On" =0, 


Ri? 2eB. 
3 —— OT, * boy, =0, 


Emin 


2 
2 E CoB = CoB b, 
Isopt Lmin RC T Onl; 2 Onl; ? 


min 


R= 2685. Вы ly +С EXE r, 
BP? B N'k o 2С N'k |? 


z(t 
Es sare ( ) Oy. 
gd( B^ / AC + A) 
Для примера провели расчет расхода топлива для автомобиля CO следующими характеристиками: полная 
масса автомобиля т = 1800 кг; коэффициенты внешней характеристики двигателя А = 0,64; Б = 1,36; С = 1,0; 


ит = 


максимальная мощность двигателя тах = 100 кВт; угловая скорость вращения коленчатого вала, 
соответствующая максимальной мощности двигателя оу = 576 c^; радиус колеса ry = 0,34 м; КПД трансмиссии 


тр = 0,85; передаточное число главной передачи 1, —3,4 и передаточные числа коробки переключения 


передач 15 = 1,0; ц = 1,5; в = 2,2; i =3,2; i = 4,8; коэффициент вращающихся масс 6 = 1,0; коэффициент 
сопротивления движению y = 0,12. 

Приведенный пример расчета расхода топлива для автомобиля с заданными характеристиками наглядно 
показывает решение задачи синтеза управления транспортным средством. Графики разгона автомобиля со 
ступенчатой коробкой переключения передач, иллюстрирующие решением для оптимального закона изменения 
передаточного числа, приведены на рис. 1. Для сравнения представлен график для разгона транспортного средства 
со ступенчатой коробкой переключения передач, который был построен по полученной аналитической зависимости. 


На первом этапе разгона до скорости X,,,, кривые совпадают. В данном случае X yy = 6.14 м/с или 22,1 км/ч. 


мин 


y, uo exp] ifi? +4ас, |-у, 
07 


– ү, 
(1) = се; | (21) 
i i? exp ЕЕ + 4а;с; | 
Vo; ш Va 
a; = Афу m G; b, _ Bi) max G. (22) 
б 
Г, м/с 
40 
30 
20 
10 
tc 
0,0 2,5 5.0 7,5 


Puc. 1. Графики разгона автомобиля co ступенчатой коробкой переключения 
передач (сухая укатанная дорога а = 100; f= 0,025; y = 0,118): 
орг — оптимальное передаточное число трансмиссии; 
icm — передаточное число, определяемое по времени разгона 
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19 дтах b; t үк +4a;c; 
cj = EN п 


i 2с. 


1 


; Q3) 


где Флтах — коэффициент мощности двигателя; Vo; — начальная скорость машины Ha і-той передаче; 
1 =1,3 — индекс, обозначающий номер передачи; V(t) — скорость машины Ha і-той передаче. 

Обсуждение и заключение. Применение авторами аналитического метода для определения передаточного 
числа трансмиссии, а также использование данного метода в практических расчетах для автомобиля с 
заданными характеристиками показывает возможность решения задачи синтеза управления транспортным 
средством с применением математического аппарата. 

Поставленная авторами работы цель — определение оптимального передаточного числа, позволяющего 
решить задачу синтеза управления транспортным средством для снижения расхода топлива — была достигнута. 
Результатом проведенных исследований являются аналитические зависимости, полученные для двух 
предельных случаев движения автомобиля. 

Проведенные авторами расчеты для первого случая показали оптимальное значение передаточного числа 
трансмиссии, которое определили в соответствии с формулами (15) и (16). Значение передаточного числа 
трансмиссии обратно пропорционально скорости движения автомобиля. Значение скорости, определяемое для 
оптимального передаточного числа трансмиссии будет меньше, чем передаточное число, определяемое по 
времени разгона. Формирование эмпирической зависимости во втором случае при переменной подаче топлива 
Kw = f(t) и оптимальное передаточное число трансмиссии были определены в результате численного решения 
уравнения (14). 

Важность проведенного исследования заключается в упрощении решения задачи синтеза управления 
неустановившимся движением транспортного средства, что наглядно показывают полученные графические 
зависимости. Подобный результат возможен при применении объединенного принципа максимума. 
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Содержание и дисперсность ферросплавов в обмазке 
при микродуговом легировании стали 
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Аннотация 

Введение. Основным недостатком традиционных процессов диффузионного поверхностного упрочнения 
стальных изделий является большая продолжительность, поэтому проблема интенсификации таких процессов 
является актуальной. Для ее решения предложено применение высокоэнергетического воздействия на материал, 
позволяющего получить упрочненный поверхностный слой из обмазки, в состав которой входят порошки 
ферросплавов, содержащие легирующие элементы. В литературе отсутствуют данные о необходимом 
содержании и дисперсности таких порошков в составе обмазки. Цель исследования — выбор размера частиц 
ферросплавов и их концентрации в обмазке для достижения максимально эффективного упрочнения 
обрабатываемого изделия. 

Материалы и методы. Для экспериментальных исследований использовали цилиндрические образцы из стали 
20 диаметром 12 MM и длиной 35 MM, на поверхность которых наносили легирующую обмазку, содержащую 
порошки ферросплавов и электропроводный гель в качестве связующего. После этого образцы погружали 
вертикально на половину длины в металлический контейнер, который далее заполняли угольным порошком с 
размером частиц 0,4—0,6 мм. Затем пропускали электрический ток величиной от 2,5 до 3,0 А в цепи источник 
питания — контейнер — угольный порошок — образец. Продолжительность процесса составляла 2—8 мин. 
Результаты исследования. Выполнена расчетная оценка электрической проводимости угольного порошка, 
рассчитаны теплофизические параметры микродугового нагрева стали: мощность, выделяемая электрическим 
током на поверхности стального изделия, плотность теплового потока, энергия единичного микродугового 
разряда. Получены выражения для расчета размера частиц порошка ферросплавов, а также экспериментальные 
зависимости толщины диффузионного слоя от размера частиц ферросплавов и их содержания в обмазке. 
Обсуждение и заключение. Результаты исследования позволили определить диапазон размеров частиц 
ферросплавов и их содержание в обмазке, позволяющие обеспечить наиболее эффективное легирование 
поверхности обрабатываемых изделий. Полученные данные будут использованы при разработке 
технологических процессов поверхностного упрочнения стальных изделий. 


Ключевые слова: химико-термическая обработка, энергия микродуги, диффузионное насыщение стали 
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Abstract 

Introduction. The main disadvantage of traditional processes of diffusion surface hardening of steel products is its long 
duration. Therefore, the problem of intensification of such processes is relevant. To solve it, the use of high-energy effects 
on the material is proposed, which allows us to obtain a hardened surface layer from a coating composed of ferroalloy 
powders containing alloying elements. There is no data in the literature on the required content and dispersion of such 
powders in the composition of the coating. The aim of this study was to select the particle size of ferroalloys and their 
concentration in the coating to achieve the most effective hardening of the processed product. 

Materials and Methods. For experimental studies, cylindrical samples made of steel 20 with a diameter of 12 mm and a 
length of 35 mm were used. On the surface of these samples, an alloying coating containing ferroalloy powders and an 
electrically conductive gel as a binder was applied. After that, the samples were immersed vertically for half their length 
into a metal container, which was then filled with carbon powder with a particle size of 0.4—0.6 mm. Then an electric 
current of 2.5 to 3.0 A was passed in the circuit power source — container — carbon powder — sample. The duration of 
the process was 2-8 minutes. 

Results. The calculated estimation of the electrical conductivity of coal powder was performed, and the thermophysical 
parameters of microarc heating of steel were calculated. These include the power released by electric current on the 
surface of the steel product, the density of the heat flux, and the energy of a single microarc discharge. The expressions 
for calculating the particle size of ferroalloy powder were obtained, as well as the experimental dependencies of the 
diffusion layer thickness on the particle size of ferroalloys and their content in the coating. 

Discussion and Conclusion. The results of this study have allowed us to determine the size range of ferroalloys and their 
content in the coating. This information is essential for optimizing the alloying process and ensuring the most efficient 
surface hardening treatment for steel products. The data collected will be used to develop improved technological 
processes for the surface hardening process, leading to improved product quality and performance. 


Keywords: chemical-thermal treatment, microarc energy, diffusion saturation of steel 
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Введение. Продолжительность различных процессов химико-термической обработки (XTO) стали может 
достигать 8—10 часов [1,2]. Поэтому ранее были предложены методы их интенсификации, основанные на 
применении высокоэнергетического воздействия на материал с использованием плазменного нагрева [3], 
электро-химико-термической обработки [4], комбинированных способов нагрева [5], лазерной обработки [6], 
микродугового оксидирования [7], нагрева в электролите [8]. Однако эти методы сложны для практического 
применения и поэтому не получили широкого распространения. Для интенсификации процессов ХТО 
целесообразно применение микродугового поверхностного легирования, которое позволяет достичь 
значительного (в десятки раз) ускорения процесса диффузионного насыщения, и не требует использования 
сложного оборудования. Особенностью данного метода является погружение обрабатываемого изделия с 
нанесенной на его поверхность обмазкой в угольный порошок с последующим пропусканием через них 
электрического тока. Интенсификация процесса насыщения достигается за счет образования ионов диффузанта 
и их последующей ускоренной диффузии в материал под действием электрического и температурного полей. В 
ранее проведенных исследованиях достигнуты положительные результаты применения данного метода для 
диффузионного насыщения хромом, молибденом, ванадием [9], вольфрамом [10]. Для комплексного насыщения 
на поверхность обрабатываемого изделия предварительно наносилась обмазка, в состав которой входил порошок 
ферросплава, содержащего диффузант. Эти исследования позволили определить структуру и фазовый состав 
полученных покрытий, но задача достижения максимальной эффективности упрочнения, для решения которой 
необходимо определение оптимальных размеров частиц ферросплава и его содержания в обмазке, в ранее 
проведенных исследованиях не ставилась. 
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Цель исследования — определение условий максимально эффективного упрочнения стальных изделий 
методом микродугового легирования за счет выбора размеров частиц порошков ферросплавов и их концентрации 
в используемой обмазке. 

Материалы и методы. Для экспериментальных исследований использовали образцы из стали 20, содержащей 
0,2 мас. % C, диаметром 12 мм, на поверхность которых наносили обмазку, приготовленную на основе 
электропроводного геля с добавкой порошков ферросплавов. Образцы погружали вертикально на половину длины в 
металлический контейнер диаметром 35 мм, который далее засыпали угольным порошком с размером частиц 
0,4—0,6 мм. Далее через образец пропускали электрический ток величиной от 2,5 до 3,0 А в течение 2—8 мин. 

Для приготовления обмазки использовали порошки феррохрома, ферромолибдена, феррованадия, 
ферровольфрама с размером частиц 0,40—0,50 мкм и связующее B объемном соотношении 1:1. 

Образцы подвергали шлифовке и полировке по стандартной методике с последующим травлением реактивом 
Ржешотарского. Микроструктуру диффузионного слоя исследовали с помощью микроскопа Neophot-21. 

Для измерения насыпной плотности угольного порошка использовали методику по ГОСТ 32558—2013. 
Порошок рассеивали по размерам частиц с помощью набора сит по ГОСТ 33 029-2014. 

Результаты исследований. Насыпная плотность угольного порошка, т.е. масса единицы объема свободно 
насыпанного порошка определена экспериментально, полученные значения приведены на рис. 1. 


П, г/смз 
0,90 T 


0,70 = = - x, мм 
0,2 0,4 0,6 0,8 1,0 


Рис. 1. Зависимость насыпной плотности П угольного порошка от размеров Х его частиц 


В диапазоне размеров от 0,4 до 0,6 мм среднее значение насыпной плотности можно принять равным 
0,81 r/cw?. Полученное значение позволило определить объемную плотность угольного порошка как отношение 
насыпной плотности к плотности каменного угля, значение которого заимствовано из справочной 
литературы [11] и принято равным 1,6 г/см}, откуда объемная плотность принята равной 0,81/1,6 = 0,51. 
Полученное значение примерно соответствует плотности упаковки простой кубической решетки (0,52). 
Поэтому, согласно [12], можно принять, что общее сопротивление порошковой среды представляет 
собой сумму сопротивлений последовательных слоев, состоящих из параллельных цепочек из 
контактного сопротивления Ак между частицами и сопротивление Ам самих частиц (рис. 2). 


Рис. 2. Электрические контакты угольных частиц в порошковой среде: 


Rk— сопротивление межчастичных контактов, Ru — сопротивление отдельной угольной частицы 


Расчет по данным рис. 2 позволил оценить общее электрическое сопротивление R порошковой среды: 


РЕ Y +Ry)d° 

= 2nH (r + dj) 

где т — количество слоев угольных частиц, r — радиус образца; d — диаметр частиц, Н — глубина погружения 
образца, Ак — сопротивление межчастичного контакта; Ам — сопротивление отдельной угольной частицы. Для 
расчетов были приняты следующие значения: d = 0,5:10° м, H = 1,5:10? м, r = 6 мм, Ак = 6 кОм, Ry = 0. Значение 
В принято равным 15,5 Ом. Полученное значение позволило рассчитать: мощность электрического тока Р = РВ, 


(1) 


Степанов М.С. и др. Содержание и дисперсность ферросплавов в обмазке при микродуговом легировании стали 


энергию единичного микродугового разряда Qo, поверхностную плотность тока f, плотность 4 теплового потока. 
Величина Qo определена как отношение мощности тока к суммарному количеству микродуг, которые в течение 
1 секунды одновременно воздействуют на нагреваемую поверхность. По данным видеозаписи значение M 


принято равным 22 500. Результаты расчетов приведены в таблице 1. 


Таблица 1 
Значения теплофизических параметров микродугового нагрева 
ТА f 102?, A/cm? P-10?, Br а: 10°, Вт/м? Qo: 102, Дж 
2,50 45 0,97 1,71 4,3 
2,75 49 1,17 2,07 52 
3,0 53 1,39 2,47 6,2 


Далее рассчитали максимальный размер алах частицы ферросплава, ДЛЯ парообразования которой достаточно 


энергии одной микродуги. 
Необходимое количество теплоты О, с учетом удельной теплоты сублимации À, можно записать как: О = Ат, 


где т — масса частицы. 


Для сферической частицы плотностью р: 


m=pV= ш а 
6 
Поэтому искомое условие имеет вид: 
Apnd? / 6€ Q). 
Максимальный диаметр частицы, необходимый для выполнения этого неравенства: 
d лах =3 z- * (2) 
пр 
Рассчитанные значения dmax представлены в таблице 2. 
Таблица 2 
Результаты расчетов 
Ферросплав FetCr Fe+Mo Fet+V Fe+w 
0,45 A/cm? 0,045 0,043 0,044 0,038 
dmax, MM 0,49 A/cm? 0,048 0,046 0,047 0,040 
0,53 A/cm? 0,050 0,048 0,050 0,042 


Таким образом, установлено, что максимальный размер частиц ферросплавов, используемых в составе 
обмазки, не должен превышать 38 мкм. 

Экспериментально определена толщина диффузионного слоя в зависимости от объемной доли частиц феррохрома 
в обмазке и ee толщины. Для этого в связующее добавляли порошок феррохрома ФХО10А с размером частиц 
40—50 мкм в количестве, необходимом для получения его объемного содержания в пределах 10-50 %. 

Результаты представлены на рис. 3. Видно, что максимальная толщина покрытия образуется при объемной 
концентрации частиц феррохрома в обмазке 50% и ее толщине 0,5 мм. Таким образом, максимальная 
интенсификация процесса диффузионного насыщения достигается при условии наибольшего содержания частиц 
ферросплавов в обмазке, а воздействию микродуг должны подвергаться частицы ферросплава, прилегающие к 
насыщаемой поверхности. 


H, мкм 


100 } 


50 — 


0,5 1,0 1,5 2,0 


Рис. 3. Зависимость глубины Н упрочненного слоя от толщины В слоя обмазки для различного содержания порошка 
феррохрома (об. %) в обмазке: 1 — 10; 2 — 30; 3 — 50 
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Обсуждение и заключение. Выполненная расчетная оценка электрической проводимости угольного 
порошка, используемого для микродугового легирования стальных изделий, позволила установить влияние 
величины электрического тока в цепи на теплофизические параметры процесса: плотность теплового потока на 
нагреваемой поверхности, мощность, выделяемую на нагреваемой поверхности, энергию единичного 
микродугового разряда, возникающего между нагреваемым изделием и окружающей порошковой средой при 
протекании электрического тока. Полученные результаты позволили определить диаметр частиц ферросплавов, 
используемых для микродугового легирования, их объемную долю в обмазке, которые обеспечивают достижение 
наибольшей интенсификации процесса микродугового легирования. Результаты исследований будут 
использованы при разработке технологических процессов поверхностного упрочнения стальных изделий 
методом микродугового поверхностного легирования. 
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Аннотация 

Введение. Изменение свойств сталей с бейнитной структурой имеет практический смысл, т. к. при бейнитном 
превращении под воздействием магнитного поля возможно улучшение пластичности стали при сохранении или 
повышении ее прочностных показателей. Научные изыскания в этой сфере касались вопросов влияния 
магнитного поля на термодинамику и смену сценария фазового превращения. Однако в открытых источниках 
нет детального описания воздействия магнитного поля на структуру и свойства продуктов промежуточного 
бейнитного превращения. Цель работы — исследование особенности влияния внешнего магнитного поля на 
сценарий и кинетику фазового превращения стали. 

Материалы и методы. Исследование проводилось на образцах из стали 65Г. Их химический состав 
контролировали при помощи оптико-эмиссионного спектрометра Magellan Q8. Термическую обработку 
(резистивный нагрев) проводили в установке для высокотемпературных исследований «ИМАШ 20-75». 
Температура нагрева — около 1000 °С, время выдержки — 10 минут. Образец охлаждали при помощи 
водоохлаждаемых электроконтактов. Внешнее магнитное поле напряженностью 400кА/м и 800 кА/м 
создавалось электромагнитом, интегрированным в вакуумную камеру установки. 

Результаты исследования. Эксперименты подтвердили возможность смены сценария превращения с 
перлитного на бейнитный при воздействии внешним магнитным полем до 1 МА/м. Получены изображения 
микроструктуры и поверхностного рельефа образцов после охлаждения в магнитном поле. Проанализированы 
кинетические изменения и зависимости объемных скоростей превращения от времени изотермической 
выдержки. Установлено, что действие постоянного магнитного поля напряженностью 1,6 МА/м увеличивает 
объемную скорость превращения в 1,808 раза (для стали 65Г) и в 1,687 раза (для стали 45Х). 

Обсуждение и заключение. Зафиксированы результаты наблюдений за изменением поверхностного рельефа 
при охлаждении без магнитного поля и в магнитных полях различной напряженности. Это позволило сделать 
вывод о стимуляции внешним магнитным полем бейнитного превращения вместо исходного перлитного. 
Микроструктурные изменения объясняются влиянием поля на магнитное состояние исходной фазы. 


Ключевые слова: магнитное стимулирование бейнитного превращения, перлитное превращение, улучшение 
свойств сталей, микроструктурные изменения сталей, вакуумное травление, магнитострикционные деформации 
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Abstract 

Introduction. It makes practical sense to change the properties of steels with a bainite structure, as with bainite 
transformation under the influence of a magnetic field, it is possible to improve the ductility of the steel while 
maintaining or even increasing its strength. Scientific research in this area has focused on the influence of the magnetic 
field on thermodynamics and on the change in the phase transformation scenario. However, there is no detailed 
description in open sources of the effect of a magnetic field on the structure and properties of the products of 
intermediate bainite transformation. The aim of the work is to study the peculiarities of the influence of an external 
magnetic field on the scenario and kinetics of phase transformation of steel. 

Materials and Methods. The study was conducted using samples made of 65G steel. Their chemical composition was 
monitored using a Magellan Q8 optical emission spectrometer. Heat treatment (resistive heating) was carried out in an 
IMASH 20-75 installation for high-temperature research. The heating temperature was approximately 1000 degrees 
1000°C, and the holding time was 10 minutes. The sample was cooled down using water-cooled electrical contacts. An 
external magnetic field with a strength of 400 kA/m and 800 kA/m was created by an electromagnet integrated into the 
vacuum chamber of the installation. 

Results. The experiments confirmed the potential for altering the transformation pathway from pearlite into bainite in the 
presence of an external magnetic field of up to 1 MA/m. Images of the microstructure and surface relief of samples after 
cooling in a magnetic field were obtained. Kinetic changes and dependencies of the volumetric transformation rates on the 
duration of isothermal exposure were analyzed. It has been found that exposure to a constant magnetic field of 1.6 MA/m 
increased the volumetric transformation rate by 1.808 times (for 65G steel) and by 1.687 times (for 45H steel). 

Discussion and Conclusion. The results of observations of changes in the surface relief during cooling in the absence 
of a magnetic field, and in magnetic fields of various strengths, were recorded. This has allowed us to draw the 
conclusion that the external magnetic field stimulates the bainitic transformation instead of the original pearlitic one. 
Microstructural changes can be explained by the influence of the magnetic field on the initial phase magnetic state. 


Keywords: bainite transformation magnetic stimulation, pearlite transformation, improvement of steel properties, 
microstructural changes in steels, vacuum etching, magnetostrictive deformations 


For citation. Dolgachev YuV, Pustovoit VN, Nefedov DV. Stimulation of the Bainite Transformation Scenario by an 
External Magnetic Field. Safety of Technogenic and Natural Systems. 2024;8(1):88—96. https://doi.org/10.23947/2541- 
9129-2024-8-1-88-96 


Введение. Как известно, магнитное поле влияет Ha термодинамику [1], в частности Ha равновесную 
температуру фазовых переходов [2]. Значительного смещения критических температур можно достичь лишь 
при использовании достаточно сильных (импульсных) магнитных полей [3], что сопряжено с определенными 
техническими трудностями при реализации в производственных условиях. Опыты свидетельствуют о 
появлении а-фазы в конструкционных и инструментальных сталях под действием полей напряженностью 
1,6-2,4 MA/M [1]. Эти процессы проходят при температурах, которые намного превышают равновесную и не 
укладываются в предполагаемый по термодинамическим оценкам сдвиг [4]. Влияние внешнего магнитного 
поля при термической обработке стали выражается в увеличении темпа [5] и дисперсности продуктов 
превращения [6]. Из [7] известно, что при промежуточном (бейнитном) превращении под влиянием магнитного 
поля напряженностью до 2 МА/м структурные изменения позволяют получать больший запас пластичности при 
сохранении или некотором повышении прочностных показателей. 

Имеющиеся в у-фазе области с ближним магнитным порядком при действии внешнего магнитного поля 
испытывают магнитострикцию. Это происходит из-за магнитной неоднородности аустенита [8] и вызывает 
изменение в поле упругих сил решетки, причем энергия образования критического зародыша @а-фазы 
уменьшается. Внешнее магнитное поле увеличивает количество и размеры ферромагнитных кластеров [1]. В 
результате при охлаждении в магнитном поле мультипликативно растет число центров зарождения. 

В [9] показано, как ближний магнитный порядок в аустените влияет на смену сценария фазового 
превращения. Из [10] известно, что магнитное состояние у-фазы определяет превращение в тот или иной 
продукт ферромагнитной о-фазы. Это может быть феррит, перлит, бейнит или мартенсит. 
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Цель представленной работы — экспериментальная проверка возможности смены перлитного сценария 
превращения на бейнитный при воздействии внешним магнитным полем до 1 МА/м. Кроме того, следует 
оценить изменение кинетики бейнитного превращения при наложении поля. 

Материалы и методы. Использовались образцы одной плавки из стали 65Г. Их химический состав 
контролировался Hà оптико-эмиссионном спектрометре Magellan Q8 (таблица 1). 

Таблица 1 


Среднее содержание элементов в образцах 
Массовая доля, % 
| C | м | c | 
65T 0,65 0,20 0,97 0,21 0,009 0,0012 0,08 0,13 


Марка стали 


Полированные образцы помещали в вакуумную камеру установки для высокотемпературных исследований 
«ИМАШ 20-75». Данное оборудование обеспечивает разрежение в рабочей камере 1,3:10%° + 6,6-10-° Па. Это 
позволяет реализовать метод вакуумного травления. Структура выявляется в результате испарения в вакууме 
при высокой температуре под влиянием поверхностного натяжения. При охлаждении в процессе фазовых 
превращений, сопровождающихся сдвиговыми процессами, на поверхности образца появляется 
соответствующий рельеф. 

Образцы подвергали резистивному нагреву до ~ 1000 °С в течение 10 минут и охлаждали с отводом тепла в 
водоохлаждаемые медные электроконтакты со скоростью ~28+32°C/c. В соответствии с диаграммой 
изотермического распада аустенита стали 65Г получаемая скорость охлаждения соответствовала пересечению 
«носа» области начала перлитного превращения. Высокая температура нагрева способствовала росту 
аустенитного зерна и вакуумному травлению его границ, что требовалось для видеозаписи поверхностного 
изменения рельефа в процессе фазового превращения. 

Для измерения температуры спай термопары приваривали в середине образца. Во время охлаждения 
протекающие на поверхности процессы регистрировали с помощью цифровой окулярной насадки Eakins. 
Полученные данные о микроструктуре поверхностного рельефа и внутри образца обрабатывали для 
определения объемной доли структурных составляющих. С этой целью задействовали аналитическое 
программное обеспечение SIAMS 800 на 16-25 полях зрения. 

Внешнее магнитное поле создавали электромагнитом, интегрированным в вакуумную камеру 
установки [11]. В экспериментах задействовали поля напряженностью 400 кА/м и 800 кА/м. Образцы без поля 
обрабатывали тем же способом, но при извлеченном из вакуумной камеры электромагните. 

Результаты исследования. При бейнитном или мартенситном сдвиговом характере превращения на 
поверхности отполированного образца должен появляться рельеф. Если этого не происходит, речь идет о 
перлитном превращении. В условиях эксперимента скорость охлаждения была недостаточна для закалки на 
мартенсит стали 65Г. Можно было ожидать конкуренции между перлитными и бейнитным превращением в 
зависимости от наличия или отсутствия внешнего магнитного поля в процессе охлаждения. 

На рис. 1 представлены скрины видеокадров изменения поверхностного рельефа, наблюдаемого при 
превращении переохлажденного аустенита в различные моменты времени (т). Каждой строке на рис. 1 
соответствует указанная напряженность внешнего магнитного поля (Н). 
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Puc. 1. Изменение поверхностного рельефа BO время охлаждения в зависимости OT напряженности внешнего магнитного 
поля (Н, кА/м) и времени превращения (т, мс) 
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На рис. 2 показаны микроструктуры поверхностного рельефа и внутри образцов после различных режимов 
обработки. 


Рис. 2. Микроструктура поверхностного рельефа (a, б, в) и внутри образца (2, 0, е) после обработки: a, г — без поля; 6, 0 — в 
магнитном поле напряженностью 400 кА/м; д, е — в магнитном поле напряженностью 800 кА/м 


При обработке без поля изменения поверхностного рельефа наблюдались лишь в отдельных мелких 
образованиях (1-я строка, рис. 1 и рис. 2 а). На видеозаписи охлаждения зерен аустенита было чуть заметно 
распространение волнового процесса по поверхности. По-видимому, это отражение фазового превращения. 
Микроструктурный анализ внутренних слоев образца показал наличие подавляющего количества (98 %) 
перлитной структуры (рис. 2 2). 

С наложением внешнего магнитного поля поверхностный рельеф формировали зарождающиеся и растущие 
бейнитные иглы (2-я строка, рис. 1 и рис. 2 6). Эти кристаллы нельзя отнести к мартенситным из-за их 
медленного роста (рис. 1). С увеличением напряженности поля повышалась интенсивность сдвигового 
превращения на поверхности (3-я строка, рис. 1 и рис. 2 в), образовывался пакетный бейнит. 

При действии магнитного поля напряженностью 400 кА/м объемная доля бейнита увеличилась 
втрое — до 6 % (рис. 2 д) по сравнению с обработкой без поля. В магнитном поле напряженностью 800 KAM 
объемная доля бейнита составляла уже 8-10 %. Это позволяет сделать вывод о стимулировании смены 
сценария превращения из перлитного в бейнитное при обработке во внешнем магнитном поле. 

Если при охлаждении исходную фазу (аустенит) сохранять как можно дольше, то усиливающийся в ней с 
понижением температуры ближний магнитный порядок (магнитная неоднородность) в своей критической 
степени приведет к атермическому мартенситному сценарию превращения. Это известно из [9]. В [1] показано, 
что наложение внешнего магнитного поля в процессе охлаждения аустенита приводит к дополнительному, 
вынужденному магнитному расслоению аустенита из-за увеличения количества, размеров и времени 
устойчивого существования ферромагнитных кластеров. В условиях текущего эксперимента искусственно 
усиленный внешним полем ближний магнитный порядок B ү-фазе способствует протеканию бейнитного 
превращения вместо перлитного, естественного для данных условий охлаждения без наложения поля. 

Под действием магнитного поля значительно меняется кинетика бейнитного превращения. Об этом 
свидетельствует отмеченное в ходе эксперимента усиление бейнитной реакции при увеличении напряженности 
магнитного поля. Кроме того, в [7] рассматривалось наложение внешнего магнитного поля напряженностью 
1,6 МА/м при изотермической выдержке различных сталей. Итоги указанной работы также подтверждают 
сделанное выше заявление о кинетике превращения. 

Для выяснения механизма влияния внешнего магнитного поля на кинетику промежуточного превращения 
учитываются: 

— специфика роста бейнитных кристаллов, который зависит от скорости отвода атомов углерода OT y/o границы; 

— возникающие при превращении структурные напряжения, обусловленные изменением удельного объема 
превращающихся фаз. 
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Данный подход обусловлен тем, что процессы релаксации у межфазной границы сильно заторможены при 
пониженных температурах [12]. В такой ситуации градиент напряжений вызывает дрейф атомов углерода 
вопреки концентрационной неоднородности. Это играет определяющую роль в скорости роста нижнего 
бейнита. Перемещение атомов С ведет к снижению его концентрации в объемах у фазы по фронту роста 
пластин бейнита. В результате возникает контролируемый концентрационным градиентом диффузионный 
поток. Он направлен навстречу растущему кристаллу и уменьшает скорость его роста. Зависимости, 
описывающие эти процессы, есть в работах Л.Н. Александрова и Б.Я. Любова [13]. 

Под действием внешнего магнитного поля возникают магнитострикционные деформации. Их упругая 
энергия может вносить определенный вклад в энергию взаимодействия диффундирующих атомов с полем 
структурных напряжений. Следует ожидать усиления дрейфа углерода и увеличения скорости роста кристаллов 
а фазы. Скорость роста по механизму дрейфа в соответствии с расчетами Л.Н. Александрова и Б.Я. Любова 
описывается зависимостью: 


y - 2PD | (1) 


Cap 
RokT|1-| 9. 
Co 


Здесь D — коэффициент диффузии углерода в аустените; Rkp — радиус критического зародыша при данном 
переохлаждении; Cup и Co — соответственно концентрация углерода в а фазе и исходном аустените. Р 
характеризует энергию взаимодействия диффундирующих атомов с полем структурных напряжений, 
вызванных дилатацией = при изменении объема в процессе превращения, и находится из соотношения: 


8 
Р = пто n 2 
9 (2) 


где г. — радиус атома углерода; © — параметр, характеризующий B твердых растворах внедрения зависимость 
постоянной решетки от концентрации атомов растворенного элемента; Ё — модуль Юнга; у — коэффициент 
Пуассона; и = (Е /2). (1 — v) — коэффициент Ламэ (модуль сдвига по направлению [100]); т = (a — b): (a * b), 
где а и Ь — размеры полуосей кристалла в форме эллипсоида вращения. 

Расчет по формуле (1) проводили для углеродистой стали при: Со= 0,7; т = 0,9; v=0,3; и = 73,5 Ina; 


—131 2 
2 500; моль) = ; = 0,01; о = 0,02. Величину Rkp определяли из соотношения 


Кр = 26 / AFo. Поверхностное натяжение на когерентной границе с = 0,2 Дж/м?; AFo при T= 600 К составляет 
315 МДж/м. Значение С.ф = 0,4 определяли как концентрацию, соответствующую температуре начала 


D = 0,0999 ev( 
© 


мартенситного превращения (600 К). Найденная таким образом скорость роста кристаллов @ фазы составляет 
~ 8:1076 см/с, т. е. является величиной того же порядка, uro и Vep ~ 10% см/с, полученная экспериментально для 
стали 65Г при видеосъемке изотермического образования рельефа. 

При включении магнитного поля измерение скорости роста дало те же результаты, т. е. экспериментально 
не обнаруживается влияние магнитострикционных напряжений на интенсификацию дрейфа. 

Допустим, упругое смещение равно истинной магнитострикции парапроцесса À = 0,5: 107* [14]. В этом случае 
при Н = 1,6 МА/м и t = 400 °С напряжения от магнитострикции составляют с), = E; = 10 Мпа. Они вносят весьма 
малый вклад в энергию взаимодействия диффундирующих атомов с полем структурных напряжений. Это 
объясняется тем, что величина магнитострикции на два три порядка меньше величины дилатации 
(0,01 + 0,07 [15]) при сдвиговом y—a переходе. По этой причине оценка вклада магнитострикции по 
формуле (1) дает исчезающе малую разницу в значениях скорости роста при обработке без поля и в магнитном поле. 

Ускорение процесса образования бейнита можно фактически оценить при помощи уравнения 
А.Н. Колмогорова [16], устанавливающего связь между степенью превращения, скоростью зарождения и 
ростом центров новой фазы. Изменение объемной скорости превращения в магнитном поле из-за уменьшения 
работы образования ферромагнитного зародыша критического размера [1] можно найти по формуле: 

Он wji 


=ехр 1+ ar А 
Vo kt Af 


(3) 


где о — объемная скорость превращения (индексы 0 и Н соответственно обозначают условия обработки без поля 
и во внешнем магнитном поле); W — энергия образования равновесного зародыша; / — намагниченность; 
Н — напряженность магнитного поля; Af — удельная «химическая» движущая сила; Af" — свободная энергия 
образования одного ферромагнитного кластера [1]. 
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Согласно классической теории Л. Кауфмана и М. Коэна [17], энергия зарождения центров новых фаз равна: 


2.7.0252 
у 95 Hu (4) 
(А+ Af") 
Здесь ДГ — работа образования ферромагнитноупорядоченных кластеров B аустените; Ө — параметр, 


учитывающий влияние энергии упругих деформаций; o — поверхностное натяжение. Из [17] известно, что 
0203 = 9,92 - 10 Дж/м. 

Вычисления проводились для температур 543 К (сталь 65D) и 628K (сталь 45Х) при напряженности 
магнитного поля Н = 1,6 МА/м. Энергия образования равновесного зародыша W определялась при удельной 
химической движущей силе Af= 150 МДж/м? с учетом напряженности поля и среднего размера ферромагнитного 
кластера ~ 1,8 нм [1], для которого получена величина Af = 0,63 МДж/м3. Расчет по формуле (3) дает для 
бейнитного превращения в стали 6ST vg / vo = 1,804, a в стали 45X он / vo = 1,665. 

Зависимости объемной скорости превращения (рис.3) получены с использованием экспериментальных 
данных [7] о степени бейнитного превращения в сталях 65Г и 45Х при различных временах изотермической 
выдержки без поля и в магнитном поле напряженностью 1,6 МА/м. 
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Рис. 3. Зависимости объемной скорости превращения (вертикальная ось, %/с) от времени изотермической выдержки: 
а — для стали 65Г при температуре 375 °С; б — для стали 45Х при температуре 320 °С. Синяя кривая отражает данные по 
обработке без поля, красная — в магнитном поле. На каждой линии отмечены экспериментальные точки 


Рис. 3 позволяет рассмотреть отношение объемной скорости превращения при обработке в поле к обработке 
без поля он/о в самом начале бейнитной реакции (экстраполированное значение B т= 1 с). Для стали 65Г 
получится 1,808, для стали 45Х — 1,687. Это близко к приведенным выше теоретическим оценкам: 1,804 для 
стали 65Г и 1,665 для стали 45Х. В процессе превращения скорости меняются. Если рассмотреть самые первые 
экспериментальные точки на графиках, то: 

— для стали 65Г при T= 40 c ou / 0 = 1,75; 

— для стали 45X при t= 28 c vu / = 1,629. 

При более длительной изотермической выдержке отношения объемных скоростей уже мало коррелируют с 
расчетными. 

Заключение. Экспериментально установлено, что, задействуя внешнее магнитное поле, можно сменить 
сценарий превращения с перлитного на бейнитное. Это объясняется усилением степени ближнего магнитного 
порядка в аустените в результате увеличения количества и размеров ферромагнитных кластеров в ү фазе под 
влиянием магнитного поля. С ростом его напряжения увеличивается магнитная гетерогенность аустенита. Как 
следствие, повышается скорость бейнитного превращения. Расчетные значения увеличения объемной скорости 
превращения при действии внешнего магнитного поля хорошо согласуются с экспериментальными данными 
для стартовых скоростей бейнитного превращения. 
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